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^ ! ABSTRACT 



I The strong variability of magnetic central engines of AGN and GRBs may result in 

highly intermittent strongly magnetized relativistic outflows. We find a new magnetic 
acceleration mechanism for such impulsive flows that can be much more effective than 
\ the acceleration of steady-state flows. This impulsive acceleration results in kinetic- 

0^' energy-dominated flows that are conducive to efficient dissipation at internal MHD 

JL ' shocks on astrophysically relevant distances from the central source. For a spherical 

^ , flow, a discrete shell ejected from the source over a time to with Lorentz factor F ^ f 

and initial magnetization ctq — Bq/Attpqc? 3> f quickly reaches a typical Lorentz 

. factor F ^ Cq^'^ and magnetization a at the distance i?o ~ ctg- At this point 

the magnetized shell of width A ~ i?o in the lab frame loses causal contact with the 

\ source and continues to accelerate by spreading significantly in its own rest frame. The 

^ ■ expansion is driven by the magnetic pressure gradient and leads to relativistic relative 

' velocities between the front and back of the shell. While the expansion is roughly 

XjT^ , symmetric in the center of momentum frame, in the lab frame most of the energy and 

0^ ' momentum remain in a region (or shell) of width A i?o at the head of the flow. This 

, acceleration proceeds as F ^ (cto-RZ-Ro)^^"^ and a ^ (7^^ {R/ Ro)^^/^ until reaching 

' a coasting radius Rc ^ Ro<^o where the kinetic energy becomes dominant: F ~ ctq 

I and a ^ 1 at Rc- Then the shell starts coasting and spreading (radially), its width 

. growing as A ~ Ro{R/Rc), causing its magnetization to drop as ti Rc/R at i? > Rc- 

Given the typical variability time-scales of AGN and GRBs, the magnetic acceleration 
in these sources is a combination of the quasi-steady-state coUimation acceleration 
close to the source and the impulsive (conical or locally quasi-spherical) acceleration 
further out. The interaction with the external medium, which can significantly affect 
the dynamics, is briefiy addressed in the discussion. 
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1 INTRODUCTION 

The first questions raised by the discovery of astrophysi- 
cal jets are how they are powered, coUimated, and acceler- 
ated. Most of them - jets from young stars. Active Galac- 
tic Nuclei (AGN), Galactic X-ray Binaries, and Gamma 
Ray Bursts (GRBs), are associated with disk accretiorQ, 
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^ The only exceptions are the jets of Pulsar Wind Nebulae as 
there are no indications of accretion disks around their pulsars. 



and this suggests that accretion disks are essential for jet 
production. The astrophysical jets seem to be highly su- 
personic as many of their features are nicely explained by 
internal shocks. In the laboratory, highly collimated super- 
sonic jets are normally produced when a high pressure (and 
temperature) gas escapes from a chamber via a finely de- 
signed nozzle. However, it seems highly unlikely that such 
refined "devices" are formed naturally in astrophysical sys- 



These jets are most likely not produced directly by the pulsars 
but ins tead form downstream of the termination shock of pulsar 
winds llLvubarskvll2002l: [ Komissarov fc Lvubarsk"vll2004D . 
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tems. They would require cold and dense gas to form the 
walls of t he chamber with a ma ssive compact object in 
the center (jSlandford fc Ree3ll974 ), but such confiKuration s 



are highly unstable (jNorman et al 



198ll : ISmith et al.lll983h 



This has lead to the idea that the coUimation of astrophys- 
ical jets may have a completely different mechanism involv- 
ing a strong magnetic field. Although this magnetic field still 
needs to be confined within a channel, the conditions on its 
geometry are less restrictive. If this field is anchored to a 
rotating object, such as an accretion disk, then it naturally 
develops an azimuthal component. The hoop stress associ- 
ated with this magnetic field component creates additional 
coUimation of the fiow within the channel. Moreover, this 
leads to a magnetic torque being applied to the rotating ob- 
ject and thus a natural way of powering outflows by tapping 
the rotational energy of the central object. 

In order to produce a relativistic flow this way, the mag- 
netic energy per particle must exceed its rest energy. Thus, 
the jet plasma must be highly rarefied. Such rarefied plasma 
is naturally produced only in the magnetospheres of black 
holes and neutron stars. Moreover, the strong magnetic field 
shields these magnetospheres and prevents them from be- 
ing contaminated by the much denser surrounding plasma. 
In contrast, young stars can eject a lot of mass from their 
surface and this seems to explain why their jets are not rel- 
ativistic. Magnetospheres of accretion disks are likely to be 
heavily mass-loaded and are not able to produce relativistic 
jets for the same reason. 

It has to be stressed that magnetic fiows must still 
be collimated externally until they become super-fast- 
magnetosonic. The magnetic hoop stress can result in self- 
coUimation of the inner core but cannot prevent sideways 
expansion of the outer sheath. However, when the fiow be- 
comes super-fast-magnetosonic, the speed of this lateral ex- 
pansion becomes smaller than the flow speed along the jet di- 
rection, and the jet remains collimated. For non-relativistic 
jets the condition of passing through the fast-magnetosonic 
surface also implies almost completed acceleration of the 
fiow (50% conversion of magnetic energy into kinetic en- 
ergy). In contrast, the relativistic jets still remain Poynting- 
fiux dominated at this point and the acceleration process 
may continue well into the super-fast-magnetosonic regime. 

The issue of the efficiency of energy conversion (from 
magnetic to kinetic form) is related to the issue of subse- 
quent energy dissipation, which is required in order to ex- 
plain the observed electromagnetic emission from both the 
jets and the structures they create when they collide with 
the external medium. Traditionally, one of the most favorite 
channels of dissipating the energy of supersonic fiows has 
been the formation of shock waves. However, in the case of 
relativistic fiows this mechanism can be much less efficient 
if the flow is Poynting-fiux dominated. First of all, it is the 
kinetic energy of the flow that is dissipate4E|, and if only a 
small fraction of the total energy is in the kinetic form then 
this already severely limits the efficiency of dissipation. Sec- 
ondly, the compression ratio and hence the fraction of kinetic 



^ This applies to fast magnetosonic shocks. At a slow magne- 
tosonic shock, the magnetic energy dissipates as well and the ki- 
netic energy can actually increase. However, slow shocks are much 
less robust and harder to generate compared to the fast ones. 



energy that dissipates also decrease with increasing magne- 
tization. Thus, in order to dissipate a signiflcant fraction of 
the available energy the flow should not only become super- 
fast-magnetosonic, but it should also become dominated by 
kinetic e nergy bef ore it is shocked (iLeism ann ct al.l l2005l : 
iMimica. Giannios fc Alovll2009l : iMimica fc Alov 201d'l . 

The magnetic acceleration of relativistic flows has been 
the subject of theoretical research for decades. The main 
focus of this research has been on the models of steady- 
state axisymmetric dissipation-free fiows (the "standard 
model"). The main reason behind this is simplicity. Only 
in this case was there a hope of building a rigorous the- 
ory. Yet, even this idealized model is rather complex, and 
solutions could be found only if an additional symme- 
try, e.g. self-simi larity, or other simpli f ying condition was 
introduced (e.g. Begelman fc Lil 1 19921 : IVlahakis fc Konigl 
I2OO3I : iBeskin fc Nokhrinall2006l '). More rec ently the problem 
was approached using numerical methods (jKomissarov et ahl 
^07,,2Q09a). 

There are a number of problems with the standard 
model, which are most severe in the case of a spherical wind. 
In this case the theory predicts an asymptotic Lorentz fac- 
tor of r ~ <^o^^, where ctq ~ Bq/Ahpq(? S> 1 is the initial 
magnetization parameter, which determines the maximum 
possible Lorentz factor corresponding to a total conversion 
of the Poynting fiux int o the bulk motion kinetic energy in a 
steady-state flow (e.g.. [Coldreich fc JuliarJll970l ). This is in 
conflict with the observations of many astrophysical sources. 
In particular, the high observed values of F in many sources 
would require an extremely large initial magnetization cro 
that would in turn imply a very high asymptotic magne- 
tization, o ~ (Tq''^ 3> 1, making it impossible to achieve 
efficient shock dissipation within the outfiow. 

A potential way to overcome this problem is by resort- 
ing to collimated outfiows. This can increase the asymptotic 
value of F and reduce that of o by up to a factor of ~ ^j^f'^'^j 
where 6'jot is the asymptotic half-opening angle of the jet. 
The coUimation has to be strong enough to preserve causal 
connectivity across the flow (in the lateral direction). The 
faster the fiow and the higher its fast-magnetosonic Mach 
number becomes, the smaller its opening angle should be. 
By the time one half of the Poynting flux is converted into ki- 
netic energy (a ~ 1), the jet half-opening angle Sjct should 
not exceed Smax = 1/F, where F ~ ctq is the jet Lorentz 
factor at that time. Observations of AGN jets do indeed 
show that Sjot < 1/F ijPushkarev et al.ll200£^ . However, for 
GRB jets with F ~ 400 (or 10^ < F < 10^ '') this constraint 
gives 6imax ==; 0.14° (or 0.018° < 6i~ax < 0.57°), which is much 
smaller compared to generally accepted values of the half- 
open i ng angle, 2° < O^.t < 30° (iFrail. Waxman fc Kulkariiil 



I2OOOI : IPanaitescu fc Kumajl20oC 

In addition, the standard theory of GRB afterglow emis- 
sion can e xplain the jet-break in their light curves o nly if 
fjctF > 1 ()Rhoadslll999l : ISari. Piran fc HalDern|[l999h . Al- 
though the Swift observations show that clea r jet breaks 
are n ot as common as we used to think (e.g., iLiang et al.l 
l2008h . this might be at least partly due to observational se- 
lection effects (Swift GRBs are dimmer on average as Swift 
is more sensitive than previous missions) , and there are still 
some clear cases for jet breaks in the Swift era. Finally, 
late time radio afterglow observations, when the flow be- 
comes sub-relativistic, provide fairly robust (no longer sus- 
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ceptible to strons relativ istic beaming) lower limits (e.g., 
lEichler fc WaxmanI [ioosh on the true energy that remains 
in the afterglow blast wave at that time, of a few to several 
times lO'^^ ergs jSerger. Kulkarni fc Frail |2004| : fPrail et al.1 
Such a large true energy, together with the inferred 
energy per solid angle in the prompt gamma-ray emission 
and in the afterglow shock at early times imply that the 
initial jet half-opening angle cannot be too small (typically 
not much less than a few degrees). 

It turns out that a transition from laterally confined 
to ballistic flow is accompanied b y a relatively short phase 
of acceleration of a differe nt kind (|Komissarov et aL I l2009bl : 
iTchekhovskov et alllioogh . Such a transition may occur in 
the coUapsar model at the stellar surface. A sudden loss of 
lateral pressure support causes a sideways expansion of the 
jet. If the jet is highly relativistic at the stellar surface the 
corresponding increase in the jet opening angle is negligi- 
ble. However, a rarefaction wave propagates into the jet and 
brings it out of lateral balance. The magnetic pressure force 
accelerates the flow in the lateral direction, which results in 
a significant increase of the jet Lorentz factor, particularly in 
the outer layers of the jet. This may alleviate the SjctT ~ 1 
problem of the magnetic model. However, as soon as the 
rarefaction crosses the jet it is well in the ballistic regime 
and the acceleration is over|f| Moreover, it does not ensure 
full conversion of electromagnetic to kinetic energy. Should, 
it happen a bit too soon and the jet remains Poynting- 
dominated. Even under the best of circumstances the re- 
sultant jet magnetization is still close to a ~ 1, which is too 



high for effective shock dissipation jLeismann et al, l2005l : 
lMimica7 Giannios fc Alovll2009l : iMimica fc Alovll201ol '). 

Given the problems with this basic case, other ideas 
have been put forward. The most radical idea is to as- 
sume that relativistic astrophysical jets do not become ki- 
netic energy dominated but remain Poynting dominated on 
all scales and that the observed emission comes not from 
shocks but from magnetic dissipation cites (Blandford 2002; 
iLvutikovl |2006| ) . In the context of the present work this 
may potentially serve as an alternative to internal shocks 
in cases where for some reason the magnetization remains 
high at large distances from the source. Others propose var- 
ious ways of increasing the efficiency of magnetic accelera- 
tion compared to the basic model, e.g., via allowing non- 
axisy mmetric instabilities and randomization of magnetic 
field (|Heinz fc BegelmanI [2OO0I ) . In fact, the magnetic dis- 
sipation may also help the transition from Poynting domi- 
nated to kinetic energy do minated states (Drenkhahn 2002; 
iDrenkhahn fc Spruitll20o3 ) . 

In this work we focus on the acceleration of an im- 
pulsive (strongly time-dependent) highly magnetized rela- 
tivistic out flow, which has rece ived relatively little atten- 
tion so far. IContopouloj (|l995l ) was first to consider the 
non-relativistic case of impulsive magnetic acceleration and 
dubbed it an "astrophysical plasma gun". The relativistic 
version presents a number of qualitatively different prop- 
erties. In § [2] we present a detailed study of a simplified 



^ This is in contrast with the highly robust mechanism of thermal 
acceleration, where for an adiabatic index of 7 = 4/3 the jet 
Lorentz factor grows linearly with the jet radius, F oc R, even in 
the baUistic regime. 



test case featuring a cold and initially highly magnetized 
((To 3> 1) one dimensional finite shell (of initial width lo) 
initially at rest (at t = 0), whose back end leans against 
a "wall" and with vacuum in front of it. The initial evolu- 
tion (§ 12.11 and Appendix is described by a self-similar 
rarefaction wave traveling toward the wall and accelerat- 
ing the Poynting-dominated flow away from the wall. At 
the end of this phase, at time to « lo/c when the rar- 
efaction wave reaches the wall, the mean Lorentz factor 
of the flow is (F) ~ f^o^^- Soon after to the shell sepa- 
rates from the wall and moves away from it (§ I2.2|l . The 
shell continues to accelerate and keeps an almost constant 
width of ~ 2/0- Using both numerical (§[2]3} and analytical 
(§[221 §[32] and Appendixes [Cl [F2} methods, we flnd that 
during the second phase the mean Lorentz factor grows as 
(F) ~ (crot/to)^^^ oc t^''^. This phase ends at time tc = toag, 
when the acceleration slows down and the shell starts coast- 
ing. At this point (F) ~ ctq and cr ~ 1. In §[3]we present crude 
but simple derivations of the main results of § [2] that allow 
us to understand the underlying physics and show that the 
results are robust - not very sensitive to the exact initial con- 
figuration. The analysis of the coasting phase (§ 13. 3p shows 
that ai t > tc the shell width increases as A ~ 2lot/tc oc t 
while its magnetization decreases as cr ~ tc/t oc , result- 
ing in a kinetic energy-dominated flow. 



In § [4] we address the apparent paradox of self- 
acceleration - how can the shell keep accelerating after it 
separates from the wall? We analyze a variation of our sim- 
ple test case in which the wall is removed when the rarefac- 
tion wave reaches it (at to). At subsequent times there are 
no external forces on the system, implying that the center 
of momentum (CM) velocity or Lorentz factor (Fcm) re- 
main constant and there is no global acceleration at t > to 
in this strict sense. Nevertheless, even though we flnd that 



Fc^ 



1/2 



remains constant, the more relevant astro- 



physical quantity is the mean value of F weighted over the 
energy in the lab frame, {r)_B, and it indeed increases as 
{F)_B ~ (aot/to)^^^ But to < t < tc- In §[5]we discuss the con- 
nection between our test case and relativistic astrophysical 
flows and study the possible implications of our impulsive 
acceleration mechanism for the dynamics of GRB and AGN 
jets. We also briefly address the interaction of the magne- 
tized flow with the external medium for GRBs. Our main 
results and conclusions are presented in § [6] 



Soon after the flrst version of our paper had appeared on 
the electronic archive (http://arxiv.org/archive/astro-ph/), 
an independent study of impulsive magnetic acceler a- 
tion was published t here as well l|Lvutikovl l2010al lbl: 
iLvutikov fc Listeij[2O10l '). indicating growing interest in this 
mechanism. Where the covered topics overlap, the results 
of both studies agree very well. As to the differences, their 
study focuses on the initial phase of fast acceleration (at 
t < to) and shock formation (when instead of pure vacuum 
the shell expands into a rarefled plasma), whereas the main 
subject of our paper is the operation of the impulsive accel- 
eration mechanism after the shell separates from the "wall" 
(at t > to). 
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2 TEST CASE: EXPANSION OF A 

MAGNETIZED SHELL INTO VACUUM 

A good way of demonstrating the basic dynamics of the 
acceleration of a highly magnetized impulsive flow is to start 
with a simple example that can be analyzed analytically or 
using simple one-dimensional simulations. To this aim we 
consider for our initial conditions a uniform shell of width lo 
with high initial magnetization, ctq = Bq/A-kpqc? 3> 1, where 
Bo is the initial magnetic field and po is the initial rest mass 
density. We choose Cartesian coordinates in which the plane 
of the shell is perpendicular to the a;-axis and the magnetic 
field is aligned with the j/-axis. The right boundary of the 
shell is at a; = and the left one is at a; = — io- To the left 
of the shell is a solid conducting wall and to the right is 
vacuum. 



2.1 Self-similar rarefaction phase 

At time t = we let the shell expand into vacuum. This 
is a well known problem that describes a simple rarefaction 
wave propagating into the shell towards the wall. The self- 
similar simple wave solution to the general case with non- 
vanishing thermal pressure is described in AppendixlX] Here 
we focus only on the cold limit (with no thermal pressure; 
the equations describing this case reduce to those of the pure 
gas case with an adiabatic index 7 = 2). 

Using units where c= 1, the local wave speed is 



A : 



(1) 



where Cms is the fast magnetosonic speed as measured in the 
fluid frame|3 In our (cold) limit 



1-l-cr 



(2) 



where 



(B/rf 



is the local magnetization parameter, while B' — B/F and p 
are the magnetic field and the rest mass density, respectively, 
as measured in the fluid rest frame. The equations of one- 
dimensional motion yield 

B^ _ B_ 
(see Appendix E)) and thus 

P 



= const 



(3) 



a = ao- 



Po 



(4) 



The backward characteristics of the simple wave (where 
the wave moves in the direction opposite to that of the flow) 
are straight lines described by 



(5) 



where — x/t is the self-similar variable. Integration of the 
self-similar flow equation gives (see Eq. [A26| for 7 = 2 or 
Eqs. [X29] and for ao = 0), 



* This is simply the Lorentz transformation of a velocity compo- 
nent parallel to the relative velocity of two inertial frames. 




Figure 1. The self-similar rarefaction wave solution at t = 1, 
using units of /q = 1, po = 1 and c = 1. The initial conditions 
are a uniform state with parameters ag = 30 and tio = at 
— 1 < X < and vacuum for a; > 0. Shown are the magnetic 
field By {top left panel), the Lorentz factor F (as measured in the 
wall frame; top right panel), the local magnetization parameter, 
cr = {B')'^/AiTp {middle left panel), the flow velocity Vx {middle 
right panel), the magnetic pressure = (S')^/87r {bottom left 
panel) and the density of total energy {solid line), magnetic en- 
ergy {dashed line), and kinetic energy {dash-dotted line) as mea- 
sured in the wall frame {bottom right panel). 



1 + V 



where 
J+ 



1 + c, 



'ms,0 



1 — Cms,0 



V 1 — Cms / 

= (\/^+ Vcro + 1)' 



(6) 



; 16(70 



where the last equality holds for cro 3> 1. This equation, in 
combination with Eqs. ((2} and Q, allows to find p — p{v) 
and then Eq. ([3]) gives B = B{v). Finally, Eq. (O allows us 
to find the dependence of all flow variables on ^. 

Figure [1] shows the self-similar solution for ao = 30, in 
units where po — lo = c — 1, at time t — 1 (when the left 
front of the rarefaction wave is about to reach the wall). One 
can see that both the left and the right fronts of the wave 
propagate at very close to the speed of light. The magnetic 
field and the total energy density distributions in the ex- 
panding shell are almost uniform (except for the thin bound- 
ary layers) . This is expected as the plasma inertia is very low 
and the electromagnetic part of the solution must be close to 
the corresponding solution of the Maxwell equations. Near 
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the right front the distributions of most flow parameters ex- 
hibit large gradients. In the plots of the Lorentz factor and 
total kinetic energy density we see narrow spikes. The max- 
imum value of the Lorentz factor can be found from Eq. ([6]) 
by setting Cms ~ 0. For ctq ^ 1 we find 



2(70 



(7) 



This is already a very high speed. However, only a very small 
fraction of the fiow energy is associated with this spike and 
the mean Lorentz factor is much lower. Figure [1] suggests 
that the mean Lorentz factor must be close to that of the 
sonic point, ^ = 0, for which Eq. ((Gj gives (for ao 3> 1) 



1/3 



(8) 



More sophisticated averaging procedures (such as the 
weighted averages over the energy or rest mass) described in 
Appendix iBl give values which are only slightly higher (see 
right panel of Fig. IBlf) and show that 



(r) 



1/3 



(9) 



is a very good estimate. 



2.2 Evolution after separation from the wall 

At the time t = to = io/cms(w — 0) ~ lo (where we still 
use units of c = 1) the left front of the rarefaction wave 
reaches the wall, and then the evolution of the shell changes. 
A secondary rarefaction wave is launched from the wall and 
propagates to the right, trying to catch up with right front 
of the original wave. However, both fronts propagate with 
speeds very close to the speed of light, and the spatial sepa- 
rations separation between them changes only very slowly - 
to the first approximation it is « 2lo. At t < to the original 
rarefaction wave does not "know" about the existence of the 
wall, and therefore behaves according to the self-similar so- 
lution for a semi-infinite shell. Att > to, however, this is true 
only ahead of the refiected rarefaction wave, at s > a;*(t) or 
5 > ^*{t) = Xt(t)/t, where x^it) is the location of the front 
of the secondary rarefaction, {xt{to) ~ —lo)- At a; > x,{t) 
the fluid continues to be accelerated by the pressure gradient 
created during the initial expansion. 

At a:: < ^^(t), however, inside the secondary wave, the 
density and pressure drop very rapidly and the fluid is de- 
celerated by the strong magnetic pressure gradient that de- 
velops just behind the head of this wave. Moreover, the total 
rest mass in this region is very small and one can describe 
the shell evolution as a separation from the wall. This is 
in contrast wi th the non-relativist ic version of this problem 
considered by IContopouloj 1 19951 ). where there is no such 
separation and the flow pressure and density peak at the 
wallQ The Lorentz factor and kinetic energy, as measured 
in the lab frame, drop strongly behind the front of this wave 
and only the part of the initial fiow that is not yet affected by 
the right rarefaction wave significantly contributes to the to- 
tal energetics (see Fig. [2|. Therefore, the "typical" or mean 

^ Thus, the relativistic dynamics of magnetized shell is even 
closer to the "plasma gun" action and also brings to mind the 
hypothetical "phasers" , all too familiar to the fans of the science- 
fiction series "Star Trek". 






Figure 2. Propagation of a highly magnetized cold shell of 
plasma. The plots describe the numerical solution at time t = 20 
for the same initial data as in Fig. [T] and use the same units. 
The top panels show the magnetic field By (top left panel), the 
Lorentz factor F, (as measured in the wall frame; top right panel), 
the local magnetization parameter, cr = (i?')^/47rp (middle left 
panel), the flow velocity Vx [middle right panel), the magnetic 
pressure pm = (B')'^ /Sn (bottom left panel), and the density of 
total energy {solid line), magnetic energy (dashed line), and ki- 
netic energy (dash- dotted line) as measured in the wall frame (bot- 
tom right panel). The front of secondary rarefaction is located at 
X ~ 18. 



(as averaged over the energy in the lab frame) Lorentz fac- 
tor of the shell should behave as the fiuid Lorentz factor 
F(^t) at the front of the right rarefaction wave (or the back 
boundary of the shell) . 

In the fluid frame the front of secondary rarefaction 
moves with the local magnetosonic speed. In the lab frame 
this corresponds to 



dx, 
IT 



«(C»)+Cm4$*) 
1-f t;(5,)Cms(C*) 



(10) 



In the ultra-relativistic accelerating regime, where « ~ 1 and 
Cms — 1 (where the latter requirement insures that there 
is still plenty of magnetic energy to drive the acceleration: 
(cr) ~ o-(^, ) ^ 1), it is more convenient to work with the 



corresponding Lorentz factors, F = (1 — ?; 
(1 — Cms)~^'^^, using the approximation 

1 . 1 



and r„ 



2F2 ' ^ ■ 

Substituting these into Eq. ((6]) and Eq. ([5]) yields 



(11) 
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i ms ~ 2r ■ 



and 



(r/r„ 



(r/r„ 



+ 1 



Combining Eq. (|12|) with Eq. (|13p we then obtain 



(12) 



(13) 



(14) 



The final step is to find ^, — (t) and substitute the result 
into Eq. (|14|) . In fact, in the ultra- relativistic regime Eq. (|10p 
yields 

da;* ^ 



/3. 



8r2(e*)rg,,(c.) 



4aor(e,) 



(15) 



When r ^ (To this can be simply approximated as /?» — 
dxt/dt « 1, which gives us 



(16) 



(see Appendix [C]) . Substituting this result in Eq. (|14ll we 
finally obtain 

1/3 



r(C.) = 



8to 



V2to 



oc t 



1/3 



(17) 



As a self-consistency check we note that since F* = (1 — 
p2y\/2 ^ ^j^gj^ /3* 1 - l/2rj and equations gS]) and 
(fTT)) imply 



09^4,\ -1/3 

^7 



r. 



which upon integration of /3, yields 
3 



2to 



1-h 



2/3' 



to 



2to 
t 



(18) 



(19) 



thus confirming the validity of Eq. (|16|l for t ^ crgto. 

Therefore, in this regime the mean Lorentz factor of the 
shell follows the law (E) oc t^^^. Moreover, for t ^ to Eq. ^ 

1/3 

(Tq in agreement with the results obtained 
2.11 Thus, we may conclude that 

(20) 



gives r(^. 
in 



This regime continues until the magnetic and kinetic en- 
ergies become comparable (and Pms — 1), which implies 
{(t) ~ 1 and (F) ~ ao- This occurs at the time 



tc — tgaQ 



(21) 



after which the shell starts coasting at a constant Lorentz 
factor (F) ~ cro (as described in 13. 3|) . 

In Appendix [C] we provide an alternative derivation of 
Eq. p7p . based on the explicit solution of the self-similar 
rarefaction wave. Furthermore, analytic expressions are de- 
rived for the rest mass M*, kinetic energy i?kin, electromag- 
netic energy Eem, and total energy (excluding rest energy) 
-E*, in the region between the head of the secondary rar- 
efaction wave and the vacuum interface: ^*{t) < ^ < = 
2[o-o(l + o-o)]^^^/(l + 2ao), as a function of ^,{t). Together 
with equation (|C8|) for t{^,) these quantities can be para- 
metrically expressed as a function of the time t, and are 




t/to 



LU 
UJ 



„2 _ „3 



10 10 10 10 10 10 10 10 10 10 10 

t/to 



















ao = 1 03 













10 10 10 10 10 10 10 10 10 

t/to 



Figure 3. Evolution of the shell, corresponding to the region 
between the front of secondary rarefaction wave (5,) and vac- 
uum interface (Xh)-, for three different values of the initial mag- 
netization: tTQ = 10^ {green lines), ao = 10^ {red lines), and 
(TO = 10* {blue lines). The top panel shows the width of this 
region, A«, and the rest mass M* within it, both normalized 
by their initial values at t = to (when the original rarefaction 
wave is secondary by the wall), as a function of t/to- The mag- 
netic flux, J Bdx, has exactly the same evolution as the total 
rest mass. The inset shows these quantities as well as the total 
energy (kinetic+magnetic) Et within this region, normalized by 
its initial value at t = to, -E*o = ^EM,0 = (''■o/2)Moc^, as a func- 
tion of t/tc (where tc = a'^to); cro = 10^, 10^, 10* are plotted 
with dotted, dashed, and solid lines, respectively. The curves for 
Mf/Mo and At/A«o a^re practically on top of each other, while 
those for Et/EfO are slightly offset, indicating a slower conver- 
gence in the limit cro c»- The middle panel shows the evolution 
of E,/EfO {thick solid lines), and its decomposition into kinetic 
{dashed lines) and electromagnetic {dashed- dotted lines) energies. 
The bottom panel shows for (tq = 10"^ the evolution of the aver- 
age values (weighted by energy -Eq. of F ((P)) and <r (((t)) 
within this region {thick lines), as well as their values at the head 
of the secondary rarefaction wave (^* ; thin lines) . 
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Figure 4. Numerical and self-similar solutions for shell's "tail" . 
The left panel shows the velocity and the right panel the magnetic 
field B. The numerical solution is represented by solid lines and 
the self-similar solution (Eq. \22\i by dashed lines. The problem 
parameters and units are the same as in Fig.[T] The time is t = 40 
and the numerical solution is shifted along the x axis so that the 
wall is now located at a; = 0. 



presented in Fig. for ao = 10^, 10^, 10* (presumably cov- 
ering the range of values most relevant for GRBs) . Similarly, 
we also derive the average values (weighted over the energy, 
according to Eq. [Blj ) of F ((F)*) and a ((f)*) within this 
region, which are shown in the bottom panel of Fig. |3] for 
CTo = 10^, along with r(^*) and (t(^*). 

One can see that up to the time t ^ tc — (TqIo the shell 
width, A», and total energy _E« remain almost unchanged, 
with A* ~ 2l(i and Et being equal to the total (excluding 
rest energy) energy in the initial solution, _E«o = -Eem,o = 
{a'o/2)Moc? . At the same time, the shell's total mass (and 
magnetic flux) slowly decrease due to the gradual advance 
of the secondary rarefaction into the shell. Kt t <^ tc, (F), 
is slightly larger than r(5*) (with the same scaling oc t^^'^), 
as expected, while (a), is slightly lower than (T(^t) even for 

t » tc. 

The shell's magnetic energy is gradually converted into 
its kinetic energy 

tc = o-Qto), as (F) ((Toi/io)^''^ in this regime. At t ~ O.OSfc, 
when these energies become comparable, while o"(^*) and 
(a), drop below unity at t/tc ~ 0.086 and 0.0037, respec- 
tively. At t > tc the shell begins to experience significant 
spreading (A./A*o « 2~''^^3{t/tcy^^ at t > tc). Its total 
mass and energy significantly decrease, indicating that the 
region between vacuum and the secondary rarefaction no 
longer represents the shell evolution. 

The numerical solution presented in Fig. [5] suggests self- 
similar evolution with characteristic linear profile for the 
flow velocity, v ~ x/t, for the region between the wall and 
the shell (we will refer to this region as the shell's tail). This 
is expected in the limit where the separation between the 
shell and the wall becomes much larger compared to lo, the 
only characteristic length scale of the problem. As shown in 
Appendix [E] such similarity solution does exist, 

1 Cr „ 1 ^22^ 



t 



B 



t l-,t2 



where ^ = x/t < 1, and d are constants. Figure [4] com- 
pares the similarity solution with the numerical solution at 
t = 40 (ao = 10=^). One can see that there is a reasonably 
good agreement between them. The first equality in Eq. (|22p 
shows that each fluid element moves with constant speed. 
This implies that the kinetic energy for any section [^1,1^2] 



of the solution is conserved. However, the magnetic energy 
of such a section decreases as oc t~^ . This indicates that the 
magnetic energy is transferred along the solution towards 
^ = 1, where this solution is no longer applicable (as ^ = 1 
implies F = 00). In order to confirm this conclusion consider 
a conserved Q that satisfies equation 



dF 



^ + ^ =0 



dQ 

dt ' dx 

Next consider a fluid element bounded by a;i = 
X2 = €,2t. The amount of Q held by this element, 

Q(ei,e2,t) = j Q{x,t)dx, 

satisfles the equation 
d 



^it and 



dt 



where 



Qi^i,^2,t)^F{^^,t)-F{^2,t) 



F = F-iQ 



is the flux of Q through the boundary moving with speed ^. 
For the energy. 



2 



Stt 



8nt^ 1 - ^2 



which represents the work per unit area and time done by 
the fluid behind ^ on the fluid ahead of ^ (the force per unit 
area is simply the magnetic pressure, / = and thus 

dW — fdx = (b^ /2)vdt). This is positive and monotonically 
increasing function of ^, which implies transport of energy 
through the tail towards the shell (C = l)j in the direction 
of motion of the flow. Clearly this is due to the work done 
by the magnetic pressure during the tail's spreading. In the 
tail's head this energy is presumably converted into the ki- 
netic energy. 

At late times after most of the magnetic energy is trans- 
formed into kinetic energy, this solution may still reasonably 
describe the tail of the flow, corresponding to ballistic mo- 
tion at F ^ (To. It implies that in the tail there is approxi- 
mately equal rest mass per decade in F: dM/dluF = Ci//?, 
dM/dlnu = Cil3 and 



M(</?) = ^ln 



1 + 

1-/3 



= Cl In [F(l + /?)] 



and equal energy per unit 4-velocity, dE/du = Ci or E{< 
u) — dE/d\nu = CiW oc u, so that most of the energy is 
carried by the fastest material. That is, deep in the tail there 
is a good part of the total rest mass but a very small fraction 
of the total energy. 



2.3 Numerical simulations 

In order to test the validity of our conclusions we have car- 
ried out numerical simulations for the evolution of a cold 
flnite shell, initially highly magnetized and at rest, as it ex- 
pands into vacuum. We numerically integrate the relativistic 
magneto-hydro-dynamic (RMHD) Eqs. (|AHA2|) in the cold 
limit, where the gas pressure is set to zero. As shown in 
Appendix [D] the equations in spherical coordinates can be 
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Figure 5. Snapshots of physical quantities from the numerical simulation of the evolution of a highly magnetized shell. Top row: 
density; middle row: magnetic field; bottom row: Lorentz factor. Each column corresponds to different times. Density and magnetic field 
are normalized by po and Bq such that B^^/iinpoc^ = uo- In the third and fourth columns, in order to follow the moving shell, the 
x-coordinate is centered on the location i^pcak of the peak of the density of the shell. 



reduced to the planar case, so it suffices to find the solution 
in the Cartesian one-dimensional geometry. 
Short term evolution: first, in order to validate our ana- 
lytic treatment in i) l2.2l we used the exact same initial condi- 
tions as in our analytic test case, namely a perfectly uniform, 
cold and highly magnetized shell at rest. At t = the shell 
occupies the region —lo < x < 0, where at a; = —lo it is 
bounded by a solid wall and in the region a:: > there is 
vacuum. The magnetic field is aligned with the y-direction. 
We have used the initial magnetization of ctq = 30. 

In agreement with our analytic analysis the shell sep- 
arates from the wall at dimensionless time t Ri to when its 
thickness in the lab frame is A « 2lo. After this time the 
solution can be described as a shell of constant thickness 
A ~ 2lo followed by a low energy tail (see Fig. [51 which 
shows the solution a.t t — 20to)- In the tail of the flow, the 
velocity grows linearly with x as predicted in the self- 
similar solution (see Appendix lE)) . 

Long term evolution: next we set out to test the long 
term evolution. We used slightly modified initial conditions: 
a shell of width lo with roughly constant density and mag- 
netic field, corresponding to a constant magnetization of 
ao = 30, whose back end touches a reflecting wall on the 
left (at X = —lo) and is tapered off to vacuum with a hy- 
perbolic tangent profile on the right over a thickness lo/W. 
That is, at t = and x > —lo we have p/po = {B/BoY ~ 
[1 -tanh(10a;//o)]/2 and cr = B'^/A-Kpc^ = B^/Anpoc^ = ao. 
We use a simple second-order accurate Harten, Lax, and van 
Leer (HLL) scheme with Runge-Kutta third-order time in- 
tegration for the numerical algorithm. The resolution is 100 
cells per lo, and the Courant number is 0.25. To follow the 
evolution of the relativistically moving shell for long times 
without enlarging the grid, we implemented a "moving win- 
dow" algorithm, where all quantities are shifted to the left 



by cAfshift cells every Atahift = 200 time steps. Thus, the 
simulation frame efi^ectively flies to the right at the speed of 
light, and the left wall becomes causally disconnected from 
the main domain. The moving window algorithm turns on 
after the shell moves away from the reflecting wall by about 
70Zo- The size of the moving window domain is lO** cells 
corresponding to lOO^o- 

Fig. [5] shows the profiles of density, magnetic field and 
Lorentz factor at several times during the simulation, while 
Fig. [6] shows the evolution of the energy-weighted aver- 
age Lorentz factor (defined in Eg. IBip . (F), with time. We 
measure time and space in units of the shell crossing time 
to = lo/c and initial shell widths lo, respectively. As ex- 
pected, the evolution has several distinct phases. First, the 
rarefaction wave propagates towards the refiecting wall, as 
seen in the first two columns of Fig. [S] The right end of the 
shell accelerates, and (F) reaches a^^^ when the rarefaction 
wave crosses the shell at t = to • At this point the shell decou- 
ples from the wall. As seen in Fig. [6] at t = to the evolution of 
(F) changes to the accelerating stage which takes it beyond 
u^/^ , increasing as t^^'^. In this regime, the shell remains 
thin (~ 2Zo, third column in Fig. [5]), leaving a low-density 
tail behind. This is the "impulsive" stage, where the right 
part of the shell accelerates at the expense of the magnetic 
"exhaust" on the left. The dotted line in Fig.|6]shows the an- 
alytical expectation during this stage, (F) = {t/toY^^ , 
for the parameters of the simulation. The agreement during 
the accelerating stage is very good. 

In the saturation (or coasting) stage, which starts 
around t/to^ao, the shell starts to spread significantly in 
the lab frame (last column in Fig. [Sjl . The evolution of (F) 
deviates from the earlier t^^^ power-law and begins to ap- 
proach the asymptotic value (F) = ctq {dash-dotted line in 
Fig. [6]), corresponding to the complete conversion of mag- 
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Figure 6. Time evolution of the energy-weighted average Lorentz 
factor of the shell, showing the rarefaction wave, magnetic accel- 
eration and saturation stages. 

netic to kinetic energy in the shellEl In the far asymptotic 
regime, the moving window of the simulation which flies at 
the speed of hght begins to outrun the shell, which moves 
with finite Lorentz factor. Thus, the last points in the evo- 
lution in Fig. [6] can be affected by the fact that a signifi- 
cant fraction of the shell material is left outside the moving 
window. However, the trend for saturation is clear. Overall, 
our simulations support very well the analytical arguments 
about the rarefaction wave, impulsive acceleration and the 
saturation (or coasting) stages of the evolution of an impul- 
sive flow. We have also experimented with larger values of 
CTo = 100, 1000 of the shell. We find that the t^^'^ evolution 
is robust and is seen in both of these cases; however, we 
did not run the simulations long enough to see the ultimate 
saturation, as the saturation time is much longer, scaling as 
cj\. We also checked that the evolution is not sensitive to 
the exact shape of the initial shell. 



3 "BACK OF THE ENVELOPE" 
DERIVATIONS 

In this section we re-derive the key results of previous sec- 
tions using crude but simple calculations which help to clar- 
ify its physics. They also show that this phenomenon is 
rather generic and not very sensitive to exact initial con- 
figuration. 

3.1 Initial acceleration 

As the shell expands and the flow develops, the electromag- 
netic part of the solution closely follows that of vacuum 



^ In this asymptotic limit, (r)^ — > 1 -|- (To/2 for the mass- 
weighted average defined in Eq. ||B3|| , while the exact asymptotic 
value for the energy-weighted average 1^)e defined in Eq. iBlll 
depends on the asymptotic distribution of dM/dT, which in turn 
depends on the exact initial conditions. It is nonetheless always 

~ (TO. 



electrodynamics. The shell (electromagnetic pulse) thick- 
ness increases by a factor of two, from Iq to li ~ 2lo and 
the magnetic held decreases by a factor of two, from Bo to 
B ~ Bo/2. At the same time, an electric field E ^ B is gen- 
erated. Since the flow is still highly magnetically dominated 
the energy conservation implies 

BqIo ~ {E^ + B^)li ~ 2B^h ■ 

On the other hand the mass conservation reads 

po^o ~ pTli , 

where p and F are the characteristic ("mean") density and 
Lorentz factor. From this we flnd that 

AnpT ^ 2 

Since the fluid frame magnetic fleld B' = B/T this gives 

1 (7 ~ . 

2 

From the MHD viewpoint the shell separates from the wall 
(i.e. looses causal contact with it) when its Lorentz factor 
just exceeds that of the fast magnetosonic speed, given by 
equation ([2]), which corresponds to a 4- velocity Ums = cr^^^. 
For (T 3> 1 this reads 

r ^ 1/2 

Combining the last two equations we find the anticipated 
results that at to when the shell separates from the wall, 

(r)~ay^ and (a) ^ a^^^ 

It is easy to see that these calculations are not sensitive to 
geometry and apply equally well to planar, spherical and 
cylindrical shells with a tangential magnetic field. 

3.2 Acceleration after the separation 

After the separation from the wall the total momentum of 
the shell no longer increases and it is mainly in electromag- 
netic form. However, the shell plasma (corresponding to the 
front of the flow) continues to be accelerated by the magnetic 
pressure gradient that has developed during the flrst phase. 
(Although, in the laboratory frame the magnetic field is al- 
most uniform the magnetic pressure is given by the strength 
of magnetic field in the comoving frame B' — B/T, which 
is non-uniform.) Similarly, the plasma at the back of the 
flow (inside the secondary rarefaction wave that develops 
and propagates into the back of the shell) in decelerated by 
the magnetic pressure gradient there. 

Magnetic flux conservation implies that Bl ~ const, 
where / is the shell width. Therefore, the electromagnetic 
energy scales as 

Eem oc B^l oc , 

and thus it decreases significantly when the I increases sig- 
nificantly, say doubles its initial value of l{to) = h ~ 2k). 
Since the part of the shell carrying most of the energy has 
a spread in the Lorentz factor of the order of Ar(t) ~ r(t) 
around its typical value, T{t), it spreads such that its width 
grows as 
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where we use units of c = 1 and the last approximate equal- 
ity holds for t ^ to (factors of order unity are dropped for 
simplicity). Now, I increases significantly at the time tc when 
the two terms on the r.h.s of the above equation become 
comparable, 

tc 



h ^ to 



Since t^ is also the time when the electromagnetic energy 
decreases significantly, we know that at tc we must have 
a{tc) ~ 1 and T{tc) ~ ctq, regardless of the value of tc, 
which we want to derive here. Thus, we find that 



tc ■ 



toT^itc 



We have already derived in the previous subsection that 
r(to) ~ (^y^ and therefore if indeed F increases as a power- 
law with time t between to and tc (which is the only viable 
option) then the power law index must be 



diogr _ iog[r(t,)/r(t„)] _ \og{ot, 



2/3x 



dXogt log(tcAo) log(cro) 3 

Thus we obtain the anticipated scaling F oc t^^^ ai to <t < 
tc- Since Tito) ~ (7^^ this implies 



T(to <t<tc 



An alternative derivation is provided in Appendix lF2l Thus, 
the scalings obtained for the test case of initially uniform 
shell are in fact rather generic. 

3.3 Coasting phase and summary of main results 

Kt t > tc the flow essentially becomes unmagnetized (i.e. 
with a low magnetization, cr < 1), its internal (magnetic) 
pressure becomes unimportant dynamically, and each fluid 
element within the shell coasts at a constant speed (ballistic 
motion). The shell coasts at a typical Lorentz factor of F ~ 
(To, where the expansion of the shell during its acceleration 
stage results in a dispersion AF ~ F in its Lorentz factor 
(i.e. that of the part of the shell carrying most of the energy) 
around this value. This causes an increase in the shell width 
in the lab frame, according to Eq. ^2'i\ . where a.t t > tc the 
second term on the r.h.s becomes dominant, resulting in 



l_ 

h 



2k> 



Cc 



Cc < 1 
Cc > 1 



(24) 



where ("c = t/tc ~ R/Rc, while tc = to<Jo and Rc ~ tc are 
the coasting time and radius, respectively. Since Eem oc 
and at t > tc -Ekin ~ E = const, then a — i5EM/£kin oc 
oc t~^ . One can summarize this result in terms of the 
lab frame time or the distance x ct of the shell from the 
wall (or source), either in terms of C^c, 



(r) 



or in terms of i^o 



A/3 

o-qCc 


(7(7^ < Cc < 1 




Cc > 1 , 






Cc' 


Cc > 1 , 


t/to ~ 


R/Ro, 



(25) 



(26) 



(r) 



■ (aoCo)'/' f < Co < , 

f^O Co > cr'o , 

■ 2/3^-1/3 1 ^ A ^ 2 
c^o Co 1 < Co < CTq , 



„2 A-1 z' ^ „2 

(^0 Co Co > CTO 



(27) 



(28) 



4 SELF-ACCELERATION: A PARADOX? 

The apparent self-acceleration of the plasma shell, which was 
described in §§[2]and[31 is rather unusual and even somewhat 
perplexing. This self-acceleration reminds of the outrageous 
tall tales of Baron Munchausen, particularly the one where 
he escapes from a swamp by pulling himself up by his own 
hair (or bootstraps). In this section we try to resolve this ap- 
parent paradox, and clarify how the shell keeps significantly 
accelerating after losing causal contact with the wall. 

At the heart of the apparent paradox lies the well known 
fact that for a closed system with no external forces the 
center of momentum (CM) velocity, /3cm, remains constant. 
This is valid not only in the Newtonian regime, but also in 
special relativity, where /3cm is the velocity of an inertial 
frame, Scm, where the total momentum of the system van- 
ishes, P' — 0, as measured simultaneously in that frameQ If 
we denote the energy and momentum as measured in S'cm 
by E' and P' = 0, then in an inertial frame S in which Scm 
moves at a velocity /3cm = /3cm i, and the total energy and 
momentum are E and P, a simple Lorentz transformation 
implies = = 0, = = 0, and 

P = P.^ rcmiPL + PcmE') = Fcm/3cm£' , 



E = Fcm(£^' + PcmPL) = FcmS' , 



/3cm = . (29) 



Since in the absence of external forces P and E remain con- 
stant, as measured in frame S, so does j3cM- 

Now, for simplicity let us consider a slight variation on 
our simple test case from § [2l where at the moment the 
original rarefaction wave reaches the wall (i.e., aX t = to 
in the lab frame, which is identified with frame S here), the 
wall is removed (and replaced by vacuum). This modification 
should not have any effect on the propagation speed and 
location of the head of the secondary rarefaction wave, /3» = 
dxt{t)/dt and ^*(t) — x,{t)/t, or on the flow ahead of it, 
at ^ > £,*{t). It would affect only the region behind the 
head of the secondary rarefaction wave. Therefore, it should 
not affect the local dynamics of the "shell" (where the shell 
refers to ^*(t) < ^ < = /3max). However, in this case, at 
t > to there is immediately no external force exerted on the 
flow, and therefore its total momentum and energy are fixed 
to their values ai t — to (for the energy this was true also 
before to since the wall was static in the lab frame): 



P{t > to) = P{to) = 



Fdt = fifo 
Sir 



Mo 



y^ao{l +ao) 



^ When viewed from this frame, it is obvious that in the absence 
of any external force the total momentum P' = remains un- 
changed, so this frame remains the CM frame, and its velocity 
PcM measured in any other inertial frame remains constant. 
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E{t > 0) = Mo 



2 + O-0 



/3c K 



^cro(l + (Jo) 
2 + O-0 



(30) 



where Mo = po/o is the total rest mass (which like P, E and 
F, is measured per unit area, given the ID planar geometry). 
In terms of the Lorentz factor, 



^c^ 



Pcm) 



-1/2 



1/2 



V4 + 3(70 \/3 



< cro , (31) 



(the last two approximations hold for cto S> 1). In Appendix 
[G]we derive the same result for /3cm and Fcm by calculating 
the total momentum in a general rest frame (simultaneously 
in that frame), and then requiring that it vanishes. 
In the CM frame, Scm, the total energy is 



E 
Fcm 



V4 + 3cro 



2 



Mo ~ ^al'^Mo 



(32) 



i.e., a factor of ~ ctq ^ 1 larger than the rest energy. 
Therefore, at late times when all of the magnetic energy is 
converted into kinetic energy, the typical Lorentz factor of 



fluid in this frame must be 



1/2 



and in particular 



E' V4 + 3(70 
= Mo ^ 2 



.1/2 



(33) 



However, since P' = this implies that comparable frac- 
tions (of the order of one half) of the rest mass would be 
moving at u' ~ ctq^^ and at u' ~ — ct^^, corresponding to 
F ~ (To and F ~ 1, respectively, in the lab frame. This pic- 
ture is supported by a direct calculations in the CM frame 
(for details see Appendix [C] and in particular the discussion 
around Fig.[G2}. 

This bares a lot of resemblance to the simple mechanical 
analogy that is described in Appendix IFll of two masses, 
m, initially moving together with a Lorentz factor F in the 
lab frame, and connected by a compressed ideal massless 
spring with potential energy E'^^^ in their initial rest frame. 
The spring is then released and fully converts its potential 
energy into kinetic energy of the two masses. In our case 
we can take m = Mo/2 and £pot — E' — Mo so that the 
final Lorentz factor of each mass is F* = E' /Mo in their 
original rest frame. Their velocities are parallel and anti- 
parallel to their original direction of motion relative to the 
lab frame, denoted by subscripts and '— ', respectively. 
If we choose F = Fcm then F, — \/4 4- 3cto/2 is slightly 
larger than Fcm resulting in F+ ~ ao and /3_ ~ —5/13. 
Alternatively, we could choose F = F* 7^ Fcm so that the 
mass at the back would be at rest in the lab frame: F_ = 
1 and F+ = r^{l + P^) ^ 2rl ^ (3/2)(to ~ ao (in this 
case E is somewhat larger than in the original case since 
we fixed E' and slightly increased F). In either case the 
mass at the front ends up with F-|- ~ cro and carries all (or 
almost all) of the momentum and kinetic energy in the lab 
frame, while the mass at the back has F_ ~ 1 and carries 
a negligible fraction of the total energy and momentum. In 
the CM frame, however, the two masses have equal energy 
and momenta of equal magnitude in opposite directions. 

Thus, the "Baron Munchausen paradox" described at 
the beginning of this section is resolved as follows. First, 
while in the lab frame the typical Lorentz factor at the time 
to when the original rarefaction wave reaches the wall is 



1/3 



into)} 

significantly higher, Fcm 



the center of momentum Lorentz factor is 



1/2 



This difference may be 



attributed to a simultaneity effect: the Lorentz factor of a 

rest frame where the total momentum vanishes as measured 

1/3 

simultaneously in the lab frame at to is indeed ~ . How- 
ever, the more physically meaningful definition of the CM 
framed requires that the momentum be calculated simulta- 
neously in that frame, and this accounts for the difference. 
Second, even though Fcm remains constant, in accord with 
our Newtonian intuition, we argue that the more astrophys- 
ically relevant quantity is (F) e - the energy weighted mean 
value of F (in the lab frame), and {F)_b does increase with 
time, approaching ~ ao at late times {t > tc). This is justi- 
fied below (and in the discussion around Fig. IG2p . 

It is by now clear that /3cm and Fcm remain constant 
at t > to, while {F)_b grows with time and approaches ~ ctq 
at late times. At such late times, t^tc = toao, when all of 
the magnetic energy is converted into kinetic energy. 



J dEp J dMVp 



dMF 



P 

E 



PcM 



(34) 



One might therefore ask, why is it more relevant to take the 
energy weighted average of F, (F)e, rather than that of /3, 
{/3) E, and then derive from it the corresponding value of F, 
(1 — (/3)|;)~^^^, which approaches Fcm at late times. The 
answer is that {F)^; is more representative of the Lorentz 
factor of the material that carries most of the energy in the 
lab frame, which is the frame where all of our observations 
are made and the external medium is at rest. This can be 
seen by using the simple mechanical analogy outlined above 
of two equal masses m = Mo/2 that end up with F_ = 1 
and F-f ~ (Jq. In this case E± = F±m, and 



(r)i 



Tj- + 1 
F+ + 1 



ao 



(1 - (/?>l) 



-1/2 



F+ + 1 



= Fcm 



r+/3+ 

F+ + 1 



1/2 



(35) 



so that using {/3)_e results in Fcm ~ (F+F-)"'^''^, which gives 
too much weight to the mass that ends up at rest (F_ = 1), 
even though it carries only a very small fraction of the energy 
in the lab frame, (F+ -I- 1)~^ ~ a^^ <^ 1. On the other hand, 
{F)_B is very close to F-|_, the Lorentz factor of the mass that 
carries almost all of the energy in the lab frame. 

The situation where part of a closed system with no 
external forces is accelerated to large positive velocities at 
the expense of another part of that system, which attains 
large negative velocities, is analogous to a rocket. If the 
rocket -ffuel start at rest with no external forces, then the 
total momentum remains zero all along. The body of the 
rocket is accelerated to positive velocities while the burnt 
fuel is thrown back with large negative velocities. That is 
why we had originally dubbed the impulsive acceleration 
of a shell the "magnetic rocket" effect. The analogy is not 
perfect, however, as rocket acceleration implies a causal con- 
nection between the body of the rocket and the exhaust. In 
the case of the magnetized shell, the decelerated material 
behind the secondary rarefaction wave is causally discon- 
nected from the forward material. In the self-similar solu- 



For example, with the former hybrid definition the velocity of 
that frame changes with time and approaches the constant veloc- 
ity of the proper CM frame only at asymptotically late times. 
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tion each fluid element is accelerated by the magnetic pres- 
sure gradient towards the asymptotic value of F ~ 2ao. The 
secondary rarefaction, however limits the duration of such 
forward acceleration. As soon as it reaches this fluid element 
its forward acceleration is terminated and replaced by the 
deceleration forced by the opposite pressure gradient behind 
this rarefaction. Thus, while the head of the shell continues 
to accelerate, the center of momentum speed for the whole 
flow remains unchanged, apart from the slow increase due 
to the wall effect. 



5 DISCUSSION 
5.1 General points 

5.1.1 Impulsive versus steady-state acceleration 

The main shortcoming of the steady-state magnetic accel- 
eration which can be successfully overcome in the impul- 
sive regime is best illustrated in the case of an uncon- 
fined spherical outflow. In the steady state regime such 
a flow accelerates effectively only up to the fast magne- 
tosonic point, where V ~ and cr ~ cf^^ . After this 

point the acceleration becomes very slo w, with F increas- 
ing only logarithmically with dista nce l|Tomimatsrj 1 19941 : 
iBeskin. Kuznetsova fc Rafikovlll998t ). resulting in Foynting- 
dominated flows on length scales of astrophysical interest. 
In contrast, the impulsive magnetic acceleration allows ef- 
fective conversion of electromagnetic energy, which leads rel- 
atively quickly to a kinetic energy-dominated flow. During 
the main phase of acceleration, after the separation from the 
wall in our test case, the magnetization parameter decreases 
with distance as a oc R~^^^ (F oc B}'^^) and then during the 
coasting phase as ct oc R"^ (F ~ constant). 

The steady-state magnetic acceleration of coUimated 
flows (jets) is generally more effective, leading to higher 
asymptotic Lorentz factors and lower magnetization com- 
pared to the case of unconfined flow. However, it still leads 
to the asymptotic values of magnetization param eter o" > 1 
iKomissarov et al.ll2009al : lLvubarskvll2009l . l2010l ). When the 
pressure distribution of confining medium is a power law 
Pext oc R'" , with Q > 2 the external confinement is in fact 
still rather insufficient. In such conditions, jets quickly de- 
velop conical streamlines and do not accelerate efficiently 
afterwards as the magnetic hoop stress, magnetic pressure, 
and electric force finely balance each other. The asymptotic 
value of the Lorentz factor is F « rmn{o\y^^9~^ , cro/2), 
where 9j is the asymptotic half-opening angle of the jet, 
and the corresponding magnetization parameter is ct ~ 

,\ 2/3fl2/3 ,N 

max( 2 ctq' 6/ ,1). 

When a < 2, the shape of a steady-state flow is 
parabolic, rj oc R°'^'^ (where rj is the cylindrical radius), 
and its Lorentz factor grows as F oc rj oc i?"/^ until reach- 
ing F « cro/2 (a ~ 1) after which the acceleration becomes 
ineffe ctive again ()Komissarov et a l. 2009a; Lyubarsky 20091. 
I2OIOI ). Additional acceleration mechanisms, such as the im- 
pulsive acceleration mechanism discussed in this paper, are 
needed to produce kinetic-energy dominated flows (As we 
have already mentioned in §[T] cr ~ 1 is still too high for effec- 
tive shock dissipation.) On the other hand, for 2 > a > 4/3 
the steady-state acceleration is faster compared to the im- 
pulsive one. However, even if magnetic acceleration initially 



occurs in a steady-state fashion and then continues in an im- 
pulsive fashion, the kinetic energy-dominated regime would 
still be reached at the same distance from the source. As 
we shall see later, such a cooperation of two mechanisms is 
natural in astrophysical context. 

A related issue is the level of variability at which the 
impulsive mechanism becomes significant. The best case sce- 
nario is when short bursts of activity are separated by rather 
long quiet periods, so that the length of almost empty space 
between shells exceeds by an order of magnitude (or more) 
the shell width. Then one can expect that the collisions 
between shells effectively occur only in the coasting phase 
where practically all of the shell energy is in the kinetic 
form (see § I3.3|l . The issue of interaction between multi- 
ple shells is best addressed numerically and this is left to a 
future work. More generally, the maximum fraction of mag- 
netic energy that can eventually be dissipated at standard 
MHD shocks in a variable flow, generated via the impul- 
sive plasma acceleration mechanism, can be estimated as 
/s = ((B^) - {|-B|)^)/(-B^). Essentially, this accounts for 
the decrease in magnetic energy during the transition to 
uniform magnetic flelcjf] This shows that in weakly variable 
flows the impulsive mechanism becomes insignificant. How- 
ever, the observations of AGN and micro-quasar jets indi- 
cate that violent bursts rather than smooth variability can 
be more characteristic of their central engines and a similar 
conclusion can be made regarding the GRB jets from the of- 
ten violent variability of their gamma-ray emission. An addi- 
tional shock dissipation may occur if within a single shell the 
magnetic field rapidly alternates, like in the striped winds 
of pulsars. In such a case, the typical gyration radius down- 
stream of the shocks caused by collisions between shells may 
exceed the stripes separation, leading to fast di ssipation of 
the alt ernating component of the magnetic field (|Lvubarskvl 
l2003bl ). 

Finally, let us address the efficiency of dissipation in the 
internal shocks. If the source activity duration is tv and the 
duration of the quiet phase between successive shell ejections 
is igap ^ then the maximum fraction of the initial mag- 
netic energy that can be converted to other forms (namely 
kinetic or internal) is /b = (1 -I- tu/tgap)^^- If we define the 
mean value of cr as the ratio of the total magnetic to non- 
magnetic energies then this impfies that (cr) > ti,/tgap when 
the different sub-shells coUideo However, internal energy is 
needed in order to power the observed variable emission in 
GRBs, AGN or micro-quasars. The fraction of the kinetic 
energy that is converted into internal energy at the internal 
shocks depends on the local value of cr at the shocks (de- 
creasing with increasing cr, especially at low Mach number 
shocks). Internal shocks between different sub-shells occur 
at _Ris ~ Rctga.p/tv Rc where the mean magnetization 
of the shell is {a{Ris)) ~ i^/tgap, i.e. close to the above 
lower limit. This suggests that the efficiency of internal en- 
ergy generation in internal shocks may significantly increase 

^ A similar issue arises in the theory of striped pulsar winds, 
where smooth fast magnetosonic waves from an oblique rotator 
eventually steepen into multiple fast shoc ks and the same es timate 
can be used to estimate their efficiency llLvubarskvll2003al l . 

And before the overall radial extent of the flow increases ap- 
preciably, so that the expression fg = {(B^) — (|-B|)^)/(B^) that 
is based on a constant total volume still holds. 
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with tgap/fu. However, it could already be quite reasonable 
even for tgap ~ tv for which /b ~ 0.5 and even though 
{a{Ris)} ~ 1, the magnetization at the head of each sub- 
shell is below average (which may improve the efficiency). 



5.1.2 Effects of geometry 

While our test case problem deals with flows with planar 
symmetry, the effects of geometry are relatively minor. It is 
easy to verify that the results of "back of the envelope" cal- 
culations of §[3] remain unchanged for flows with cylindrical 
and spherical geometry. Appendix iDl shows the mathemati- 
cal reason for this - a suitable variable substitution reduces 
the equations with spherical symmetry to those with the 
planar symmetry. From the physical point of view this rel- 
ative lack of sensitivity to geometry is based on the fact 
that the key factor in the impulsive acceleration is the flow 
expansion in the direction parallel to the direction of mo- 
tion, whereas the symmetry of the flow mainly regulates the 
rate of expansion in the transverse direction. Due to the 
transverse expansion of jets the transverse magnetic field, 
which we assume to be dominating, decreases as B± cc 
whereas the specific volume increases as V (x r^, where rj is 
the transverse length scale. The specific electromagnetic en- 
ergy remains unchanged, -Eem oc B^V oc r°, and hence the 
transverse expansion does not lead to magnetic acceleration. 



5.1.3 Test case and astrophysical fl.ows 

The initial configuration of our test case problem can be rel- 
evant for eruptive astrophysical phenomena involving fast 
magnetic reconnection and restructurin g of magnetic field 
configuration, like the magnetar bursts (|Lvutikovll26ol ). In 
many other cases, an astrophysical central engine may op- 
erate rather steadily on relatively long time scales. These 
scales have to be compared with the time scale required for 
the flow, which is powered by the central engine, to reach the 
fast magnetosonic point of the steady-state solution. Once 
the jet propagates beyond this point, its inner part becomes 
much less effected by the waves which are generated at the 
jet head. In particular, if the jet expands into a relatively 
empty channel, then the rarefaction wave, which propagates 
in the comoving jet frame only with the fast magnetosonic 
speed, will be conflned to the jet head and unable to prop- 
agate upstream. Therefore the fraction of the jet affected 
by this wave will rapidly decrease in time. For such cases, 
a shell moving with super-fast-magnetosonic speed will be 
a more suitable initial configuration compared to the static 
shell next to a wall of our test problem. 

After the shell of our test problem had separated from 
the wall, the plasma acceleration was driven by the mag- 
netic pressure gradient, that had been developed in the shell 
prior to its separation. Thus, it is reasonable to investigate 
whether a similar pressure distribution can develop in the 
case where there is no wall but the shell is initially moving 
with a super-fast-magnetosonic speed. In this case two rar- 
efaction waves will be moving into the shell, one from its 
head and another from its tail. However, due to the proper- 
ties of relativistic velocity addition the head rarefaction will 
be moving across the shell much faster compared to the tail 



rarefaction which will almost "freeze" at the shell tail. In- 
deed, in the laboratory frame the tail rarefaction propagates 
with the speed 



/3t 



(36) 



where the last equality holds for F, Fms 3> 1. The length of 
the jet affected by this wave grows at the rate 

1 

2f2 ■ 



A/3t = /3t - /? : 



(37) 



The head rarefaction propagates with the speed 



/3h 



1 - PPu 



F2 + Tl 



(38) 



and the length of the jet affected by this wave grows at the 
rate 

A/3h = /3 - /3h ==; TTTT^ » . (39) 



A r n 

Thus, the head rarefaction crosses the shell first and creates 
the magnetic pressure gradient which accelerates the shell 
in the direction of the head, just like in our test case after 
separation from the wall. 

Moreover, in the rest frame of the shell the head rarefac- 
tion starts propagating much earlier than the tail rarefaction 
(due to simultaneity effects, since they start more or less si- 
multaneously in the lab frame), and this is the reason why 
the tail rarefaction covers only a very small fraction of the 
shell even though in this frame the two rarefaction waves 
propagate at the same speed (/3ms) in opposite directions. 
In this frame, at the time the two rarefaction waves meet 
(very close to the back end of the fiow), the configuration is 
very close to that of our modified test case, where the wall 
is removed when the head rarefaction wave reaches it. This 
fact is used later on, in the derivation leading to Eqs. H56|) 
and ((57)) . 

Another, important issue is whether the space between 
different ejecta in astrophysical jets can be considered as 
empty. Indeed, if the jet production is completely inter- 
rupted from time to time then the external gas may rush 
into the gaps between different ejecta. The speed of such a 
lateral fiow is obviously limited by the speed of light, and 
this sets the lower limit on the length scale at which the 
gaps can be considered as empty. 



Ctyt 



where ty is the variability time-scale of the central engine 
and 6j is the jet half-opening angle. (At this distance from 
the central source the shell's cylindrical radius is comparable 
to the length of the gaps between shells.) For GRB jets with 
t„ ^ 4 ms and 9j ~ 0.1 this gives us i?min ^ 1.2 x 10^ cm, 
and for AGN jets with t„ ~ 10 days and 9j ~ 1° this gives 



us Ru 



1.5 X 10 cm. The sound speed of the surround- 



ing gas at such distances can be much lower than the speed 
of light and one may expec t the empty gaps to app ear at 
smaller distances than i?min l|Lvutikov fc Listerll201o[ ) . How- 
ever, the ejecta will most certainly drive shock waves into the 
surrounding gas, heating it to higher temperatures near the 
jet channel. On the other hand, the increased buoyancy of 
this gas will result in an outfiow, which may become a super- 
sonic wind. This will effectively reduce the speed with which 
this gas expands into the jet openings. In fact, if 6j > l/M^, 
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where is the wind Mach number, then the wind gas wiU 
be unable to reach the jet axis. 

Moreover, the jet may have this kind of protection from 
the beginning if the accretion disk produces its own super- 
sonic wind (we assume here that the relativistic jet is driven 
by a Kerr black hole). Close to the source, where Aluj < 1, 
the disk wind may still try to fill the polar region, thus cre- 
ating an obstruction for the re-born jet. However, it could 
be pushed aside by the jet on the time scale required for 
the jet to overtake the wind, ~ t^v^jic — where 
is the wind speed. Using the cited variability scales and 
~ 0.1c we then find that empty gaps may appear al- 
ready beginning from the distance of ~ 10^ cm for GRB 
jets and ~ 3 x 10^^ cm for AGN jets. 

An impulsive operation of the central engine may well 
results in trapping of some amount of weakly magnetized 
and dense external gas in the gaps of intermittent highly 
magnetized jet. This gas will then be accelerated by the 
jet, leading to development of Rayleigh- Taylor instability, 
turbulence and mixing. Clearly, this important issue requires 
further investigation. 

5.2 Application to GRB jets 

We start by considering the propagation of a single shell 
produced during an active phase of central engine of dura- 
tion t„. This time may correspond to the whole duration of 
gamma ray burst or to the duration of one of many shells 
produced during the active phase of its central engine. The 
exact nature of the jet variability is not known. In the col- 
lapsar model for long GRBs and the binary merger models 
for short GRBs, this may be related to advection of mag- 
netic field with different polarity onto the black hole, sim- 
ilar to what has been see n in recent numerical simulations 
iBarkov fc Baushevl2009l 'l. In any case, the shortest variabil- 
ity time scale in GRBs is probably given by the "viscous" 
time of the inner disk. For neutrino cooled disks this is 




where M is th e black hole mass and a is the parameter of the 
a-disk model (iPopham et al.lll999l ). In the alternative model 
of GRB central engine, which replace a super-accreting black 
hole with a millisecond magnetar, the nature of variability 
has to be different. It could be driven by a violent restructur- 
ing of magnetar magnetosphere, e.g. rising buoyant magnetic 
loops and magnetic reconnection. A relatively mild case of 
such restructuring, with the characteristic timescale of order 
of ~ 20 ms, has been seen in recent numerical simulation s 
of magnetar driven GRB jets l|Komissarov fc Barkov|[2007^ . 
This time scale gives us one of the characteristic length scale 
of this problem, the shell width 

Z = ct„ = 3 X 10** f — cm . 

V 10 ms/ 

(We use the name shell rather loosely here to describe the 
ejecta, which can be rather elongated and better described 
as a jet close to the central engine.) There are many other 
important scales in this problem. 

As we have already commented on, the initial accelera- 
tion of the flow can proceed in a steady state fashion. This 
brings into consideration the radius of the light cylinder, 



He, the distance to the fast magnetosonic surface, -Rms, the 
distance up to which the steady-state acceleration mecha- 
nism remains effective, -Rs, the distance at which the im- 
pulsive acceleration mechanism kicks in, -Rcr.t, the coast- 
ing distance Kc, and finally the distance where the shell 
begins to decelerate due to the interaction with the inter- 
stellar medium or stellar wind gas, -Rdoc. There are many 
unknowns in this problem. In particular, it is difficult to say 
what is the exact nature of the collimating medium. The 
jet is unlikely to be in direct contact with the collapsing 
star. The hot jet cocoon and the wind from accretion disk 
are more likely candidates. Let us suppose that the external 
pressure scales as pcxt oc BT^ ^ the most favourable case for 
the steady-state collimation acceleration mechanism. Then 
the steady-state jet is parabolic, R cc (where r is the 
cylindrical and R is the spherical radius), and beyond the 
light cylinder, ric — c/Q, it s Lorentz factor increase s as 
r ~ (r/ric) ~ (R/ric)'^^^ (e.g. iKomissarov et al.ir2009al l. At 

1/3 

the fast-magnetosonic surface F « Ug and thus this surface 
is located at 

Rn,s ~ ricaf/^ . (41) 

If the jet is powered by a rapidly rotating black hole (a = 
0.9) then ric ~ iRg, where Rg = GM/c? is the gravitational 
radius of the black hole. For the typical parameters of GRBs 
this gives us 



and 




Thus, for the time scale of the central source variability 

the source will be able to produce a steady-state super- 
fast-magnetosonic fiow. Since ~ 2 ms is about the shortest 
timescale for the variability of the central engine (see the end 
of § I5.1.3|l this must be always the case and the effects of 
steady-state collimation acceleration have to be taken into 
account. 

As the jet propagates into an almost empty channel 
cleared by the previous ejections, there will be a rarefaction 
wave in its heads, making its way into the jet. However, it 
will occupying only a small fraction of the jet length. In- 
deed, in the source frame the speed of this rarefaction is 
given by Eq. pS)) . which for F^ 2> F^^ imphes /3h — 1, and 
thus the jet length grows much faster than the width of the 
rarefaction wave in its head. 

The collimation acceleration becomes ineffective when 
the jet half-opening angle, 6j, exceeds the Mach angle, ^m, 
which is given by 

r> Fms/^ms /^tr\ 

where A^ms is the relativistic fast-magnetosonic Mach num- 
ber. For F ^ Fms 2> 1 and thus ^ 1 this reduces to 
~ Fms/F so that the critical half-opening angle is given 

by 
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(46) 



where r„ 



.1/2 



CTq r 1/^. At this point the Lorentz 



factor and magnetization parameter of the jet are 

1/3 / fl, \ -2/3 



^l/3fl-2/3 



"3 

0.1 



and 



v2/3 



2/3 



(47) 



(48) 



respectively. (The analysis of flows coUimated by an ex- 
ternal med i um w ith a power-law pressure distribution by 
iLvubarskvl l|2009t ) leads to a result that differs from this one 
only by a factor of order unity.) In principle, both ao and 9j 
can be estimated from observations of GRBs and their af- 
terglows. In particular, ctq can be determined using the mea- 
surements of Lorentz factor via F < ctq, where the equality 
corresponds to full conversion of the electromagnetic energy 
into the bulk motion kinetic energy. The actual location of 
the point where the jet Mach angle reaches its critical value 
and the jet enters the freely expanding regime is less cer- 
tain as it depends on the exact pressure distribution of the 
confining medium. For pext oc we have 



Rs ~ ricFs 



2/3/1-4/3 



-4/3 



For the parameters typical for GRBs this gives us 



Rs 



1.3 X 10-' 



(49) 



(50) 



This is significantly lower compared to the radius of long 
GRB progenitors, which is believed to be of the order on the 
Solar radius, Rq « 7 x 10^" cm. Beyond Rs the coUimation 
acceleration is no longer effective. 

In order to find the scale at which the impulsive accel- 
eration comes into play we first need to estimate how long it 
takes for the head rarefaction to cross the shell. The length 
of the section of the shell which is affected by the rarefac- 
tion grows with time at the rate A/3h given by Eq. (|39p . 
corresponding to a crossing time 



vrn,J 2 Vrn,J 



(51) 



where the last equality assumes a super-fast-magnetosonic 
regime. By this time the shell will propagate the distance 

p2 



Ctv- 



For the typical GRB parameters this gives us 



Ctyb 



R. 



(52) 



(53) 



In the frame moving at a Lorentz factor Fs, at the time 
the two rarefaction waves meet (very close to the back end 
of the shell), the configuration is very close to that of our 
modified test case at the time when the head rarefaction 
wave reaches the wall and the wall is removed. The main 
difference is that the initial magnetization parameter is CTs 
and the initial shell width is TsCtv Thus, after the passage 
of the head rarefaction wave the typical shell Lorentz factor 
in this frame is 



r. 



.1/3 



{ao6 



\2/9 



(54) 



and the typical value of the magnetization parameter is 
o"cr,h ~ o's^''* ~ {ao9j)^''^. Therefore, the shell Lorentz factor 
in the lab frame is 



T.Ts 



5/9a-4/9 



™(S^)""(S 



-4/9 



.(55) 



This Lorentz factor is only one order of magnitude below 
the maximum value given by ctq. In fact, this may still be 
only a conservative estimate as we have not taken into ac- 
count the acceleration related to the transverse expa nsion of 
the jet when it crosses the stellar surface (|Komissaro v et all 
l2009al : iTciiekhovskov et al.ll2009h . This additional accelera- 
tion may well increase the mean Lorentz factor by a factor of 
few (in what follows we denote this factor of k). This brings 
the Lorentz factor up to Fcr.t ~ KFcr.h and the magnetiza- 
tion parameter down to (Tcr.t = crcr,h/«:. However, even after 
this the jet magnetization is still too high for effective shock 
dissipation. 

For simplicity we assume that the stellar radius R, 
where F increases by a factor of k is -Rcr.h < R* <^ Rcr,t, 
where -Rcr.t is the radius where the tail rarefaction crosses 
half of the original shell. We now use the similarity between 
our modified test case at io and the shell in its rest frame 
prior to the crossing of the head rarefaction (referred to as 
the shell's "initial" rest frame) at the time when the head 
and tail rarefaction waves meet, as discussed below Eq. (|39|) . 
Here the shell's "initial" rest frame would be moving with 
a Lorentz factor F/ ~ k^/^F^ relative to the lab frame, 
rather than Fs. This can be understood from the fact that 
we require that after the passage of the head rarefaction 



0"cr,h 

.3/2/ 



-1 2/3 , 

K (Js SO that af 



2/3 

F/ ~ ao/af ~ K='/2(o-o/o-,) ~ k^/^F, 

Now, recall that in our test case significant additional 
acceleration after the shell separates from the wall (i.e. at t > 
to) starts only when the head of the secondary rarefaction 
wave - identified here with the tail rarefaction wave after it 

1/3 

meets the head rarefaction - reaches fluid with F ~ CTq , 
i.e. F(^») ~ f^o^^ ■ This corresponds to ^» ~ i.e. .^J, ~ 
in the frame moving at Lorentz factor F/ (referred to as 
the comoving frame), corresponding to the middle of the 
shell after the passage of the head rarefaction wave. In our 
original test case this corresponded to a single dynamical 
time (« to, i.e. between t = to and t « 2to), so we did not 
pay attention to this. The comoving shell width at the time 
when the two rarefaction waves meet is 2TfCtv and therefore 
in our present case it takes the tail rarefaction wave a time 
t'cT,t ~ F/fi, to reach ^J, ~ in the comoving frame, which 
corresponds to a time fcr.t ~ T^ti, in the lab frame. This 
result can also be obtained using Eq. p7[l with F — F/, 
tci,t ~ tv/Af5t ~ F^t„ — >■ F^t„. This corresponds to a radius 



D 1^2 , 3 2/3/)-4/3 , 



, (56) 

2/» / e, \ -4/3 / u 



Vio3/ Vo.iy Vioms/ 

where the impulsive acceleration with F oc R^''^ begins 



cm 



(57) 



for _Rcr,t < -R < -Rc. As a consistency check we verify that 
this gives Rc ~ Rcr,\{(^o/T„,t)^ ~ i?cr,to-cr,t ~ crlctv, as it 
should from the general considerations outlined in § 13.21 
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We find that i?cr,t 3> Rs, which suggests that the 
steady-state colhmation acceleration and the impulsive ac- 
celeration are scale separated. At i? = -Rcr,t the ratio of the 
shell's cylindrical radius, Tj = 9jR, to its width, Ij = ctv, is 



3 2/3 



-1/3 



(58) 



Thus, "shell" is indeed a suitable name for the flow at the 
stage of impulsive acceleration. The coasting radius is given 
by For the typical parameters of GRBs this gives us 



Rc 



2 . 



; 3 X 10' 



( 



— 17- 



iv 



(59) 



VlOV VlOms/ 

and at R> Rc the shell coasts at F ~ cro while its magneti- 
zation rapidly decreases as a ~ Rc/R- 

Within this model, prompt gamma-ray emission due to 
dissipation in internal shocks between different shells within 
the highly variable outflow naturally occurs in the region 
R ~ (l-lO)-Rc. On the one hand, the mean plasma magneti- 
zation is (T ~ 1 at ii = 7?c and then it decreases linearly with 
the distance. Thus, one of the conditions for effective shock 
dissipation, o <^ 1, is satisfied in this region. On the other 
hand, the width of individual shells begins to grow linearly 
with the distance, allowing their collisions. (For R < Rc the 
shells keep an almost constant width.) Moreover, the vari- 
ation of the flow Lorentz factor in the coasting regime is 
rather large, AF ~ F, which may potentially help increase 
the e fiiciency of shock dissipatio n up to ^ 10% (|Beloborodovl 
I2OOOI ) or even higher faobayash ifc SarilliooH ). 

In order to test the viability of our impulsive magnetic 
acceleration mechanism, the coasting radius, Rc, has to be 
compared with the deceleration radius, -Rdoc, at which most 
of the energy is transferred to the swept-up shocked external 
medium. In the "thin" shell regime (see belo w) where Rc < 
Rdcc and F(i?doc) ~ o-o, -Rdcc is given by (e.g.. lGrano^l2005h . 



Rd 



(3 - k)Ei, 



2.5 X 10^ 



l/(3-fe) 



(60) 



-1/3^,1/3 
'0 -'^iso,53'^ 



-2/3 







1.8 X io^^A::^£is 



k 



'ISO, ssctq 3 cm 

for a spherical external rest mass density proflle poxt = 
AR-'', where ^0,3 = (^o/W^, Sj.o = W^'^Ei.^^sa is the 
isotropic equivalent energy in the flow, n = no cm~^ is 
the external number density for a uniform external medium 
(fc = 0) and A = 5 X 10^^ A, g cm~^ for a stellar wind en- 
vironment (fc = 2). In some GRBs the afterglow onset time 
is observed, which is identified with the observed decelera- 
tion time, fdec ~ ^dGc/2cF^(i?dcc), and may be used to infer 
the values of F(_Rdec) and i?dcc, typically gi ving values of 
r(_Rdec) of a few hundre d and i?dcc ~ 10^^ cm (|Sari fc PiranI 
ll999l : lLiang et al.ll2010l ). 



^dcc — 



(3 - k)Eisotdc 



2tvAc{1 - 
1.0 X lO^^r 



l/(4-fc) 



(61) 



-l/4pl/4 ,1/4 
-^iso,53''doc,2 



1.8 X 10-^-/^<.3e,, cm 



fc = 



fc = 2 



where tdcc/ {1 + z) — 100tdcc,2 s. Note, however, that this 
method has an observational bias towards low values of 
r(-Rdoc) and large values of -Rdcc that correspond to large tdcc 



values, since small tdcc values are hard to measure as optical 
or X-ray follow-up observations usually start at least tens of 
seconds after the start of the prompt gamma-ray emission. 
Nevertheless, even though F(7?dec) is the Lorentz factor of 
the shocked external medium and it is close to that of the 
original ejecta only for a Newtonian or mildly relativistic 
reverse shock (the "think shell" case, where tdcc > Tgrb, 
Tgrb being the observed duration of the gamma-ray emis- 
sion from the GRB), even for the "thick shell" case (where 
tdcc ~ Tgrb) this method gives r(i?dcc) ~ Fcr and i?dcc ~ 
Rci, which is the correct deceleration radius in this regime 
(the critical values of the Lorentz factor, Fcr, and radius, 
RcT, are provided below). In both regimes tdcc ^ ^grb , so 
using Tgrb instead of tdcc in Eq. (|6Hl gives Rct, which is a 
lower limit on the value of -Rdcc. 

Only when Rc ^ 0.1-Rdcc the internal shock mechanism 
can be sufficiently effective to explain the prompt gamma- 
ray emission. Equations (|61|l , (|59p and (I40p show that this is 
satisfied only when the characteristic variability time scale 
of the central engine is not much longer than the viscous 
time scale of the inner disk (t„ ~ lO"'^ s), even though 



for -Rdc 



10^'' cm, -Rc < 10^'' cm requires t„ < 0.3 s or 



an observed variability time 1 s for a typical redshift 
of 2 ~ 2. For long GRBs, with the mean duration in the 
source frame of about ~ 10 s, this implies between a few 
tens to about one thousand of individual shells. For short 
GRBs, with the mean duration in the source frame of about 
~ 0.3 s, this number can be reduced down to between about 
a few to a few tens. Moreover, it is generally easier to obtain 
Rc < Rdcc for a uniform external medium than for a stellar 
wind environment, since -Rdcc is typically much smaller for 
a stellar wind. 

Now we briefiy discuss the interaction with the exter- 
nal medium, and when it strongly affects the flow (a more 
detailed analysis will be presented in a separate work). For 
simplicity, we shall consider a single shell and discard fac- 
tors of order unity. Let us consider a spherical outflow of 
duration to, radial width -Ro « cto, energy E, and lumi- 
nosity L « E/to, propagating into a spherical external 
rest mass density proflle pcxt = AR~'' (with fc < 10/3). 
The regime where the -Rc < -Rcr < -Rdcc (i.e., where at 
Rc only a small fraction of the total energy is transferred 
to the shocked swept-up external medium), corresponds to 
the well-known "thin shell" (or initially Newtonian reverse 
shock) case f or the deceleration of a coasti ng unmagnetized 
{a < 1) sheU (|Sari fc Piranlll995l : ISarilll997l ). which has been 
investigated in the context of GRBs. Due to the spreading 
of the shell at -R > -Rc (because of a spread AF ~ F in its 
Lorentz factor), the reverse shock gradually strengthens and 
becomes mildly relativistic at -Rdcc (which in this regime is 
given by Eq. [SD]) where it finishes crossing the shell, and 
may produce a bright emission that peaks at an observed 
time tdcc ~ (H-z)-Rdcc/cCTo ~ (-Rdcc/-Rc)T'GRB > Tgrb (i.e., 
after the the end of the prompt GRB emission). 

For the other ordering of the critical radii, -Rdcc ~ Ra < 
Rc (which may occur for large values of to or a stellar wind 
external medium), the outflow generally never reaches a 
coasting phase (so Rc loses its physical meaning as a coasting 
radius), since the magnetized shell starts being signiflcantly 
affected by the swept-up shocked external medium when the 
latter still has only a small fraction of the total energy. The 
impulsive acceleration, F ~ {aoR/ Ro)^^^ , proceeds from -Ro 
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up to a radius Ru ~ -Rcr(o-o/rcr) "'''^^ where Fa ~ 
{ERo/Ac'y^^^^-''^ and ~ R^Vl ~ (Si^o/Ac^)!/**-*), 
Then, at Ru < R < Ra the typical Lorentz factor of the 
magnetized shell becomes similar to that of the swept-up ex- 
ternal medium, F ~ (L/Ac^)^'''*7?''="^'/*, and is determined 
by the pressure balance at the contact discontinuity that 
separates these two regions. This is a phase of either a mod- 
est deceleration (for < 2) or a reduced acceleration (for 
2 < k < 10/3) that occurs while the outflow is still highly 
magnetized (cr ^ 1 for i?dcc ^ Rc)- Therefore, there might 
not be a reverse shock going into the original magnetized 
outflow, and even if such a shock develops, then it would be 
very weak and could dissipate only a very small fraction of 
the total energy. Finally, at Rct where F ~ Fcr, most of the 
energy is transferred to the shocked external medium (so 
that in this regime - Rdec = Ra)- Therefo r e, at R > _Rcr the 
flow approaches the lBlandford fc McKed (|l976t ) self-similar 
solution for a spherical constant energy relativistic blast- 
wave going into an unmagnetized external mediuml^ At 
-Rcr the magnetization is still high, a{Rcr) ~ co/Fcr, where 
this generalized "thick shell" regime (iidcc ~ Rcr < Rc) 
corresponds to ctq > Fcr, while the "thick shell" regime 
(Rc < Rci < Rdcc) corresponds to ao < Fcr- 

Altogether, the acceleration of an initially highly mag- 
netized {ao S> 1) impulsive outflow via the impulsive ef- 
fect and its deceleration due to the interaction with the 
external medium are tightly coupled and cannot be fully 
treated in isolation. That is, the magnetic acceleration nat- 
urally sets the initial conditions for the interaction with 
the external medium, and realistically one cannot simply 
assume any arbitrary initial configuration of the magne- 
tized outflow near the deceleration radius. Moreover, in the 
highly magnetized "thick shell" regime there is an interme- 
diate phase (Ru < R < Rcr ~ Rdcc), whcrc the magnetic 
acceleration and the deceleration because of the external 
medium balance each other, resulting either in a reduced 
acceleration or in a relatively modest deceleration (as out- 
lined above). If the outflow starts highly magnetized, then 
it can decelerate either in the unmagnetized "thin shell" 
regime (with (7(_Rdcc) ~ Rc/Rdcc < 1 for ao < Fcr) or in 
the highly magnetized analog of the "thick shell" regime 
(with a{Rdcc) ~ ao/TcT > 1 for ao > Fcr). There is no high 
magnetization "thin shell" regime, and in order to be in the 
low magnetization "thick shell" regime, a single-shell flow 
cannot start highly magnetized {ao 3> 1)EEI 

5.3 Application to AGN jets 

We can apply the results obtained in the previous section to 
AGN jets simply via appropriate rescaling. First, the char- 

Tile previous regime, Ru < R < Rcr corresponds to anotiier 
variant of tliat solution, witli energy injection into the siiocked 
external medium by a relativistic wind from tlie central source. 

Tiiis might still be possible, under favourable conditions, for 
a highly variable flow with a large number of high contrast sub- 
shells that accelerate independently, and then quickly collide and 
merge into a wider shell via internal shocks, soon after reaching 
their coasting radius and cr < 1. This would amount to quasi- 
continuous energy injection by subsequent sub-shells just after 
the deceleration radius of the first sub-shell, thus increasing the 
deceleration radius of the whole flow. 



acteristic masses of black holes are hig her, ~ fO^-lO^M©. 
They are radiation cooled and the corresponding shortest 
variability time scale is 

^^,nr.^l0(ii^)"(^)days, (62) 

where 5 = Hd/Rd is the ra tio of the disk height to its radius 
dShakura fc Sunvaev|[l973l ). 

Second, the Lorentz factors of AGN jets can be mea- 
sured directly via observation of proper motion of their 
knots. Such observations (mainly VLBI radio observations) 
indicate relatively low Lorentz factors, of the order of a few 
for weak radio sources (FR-I type) an d (F) ~ 10 with a 
tail extending up to F ~ 50 for blazars l|Lister et al.ll2009l ). 
The rapid variability of gamma-ray emission from some 
AGN sug gests the possibility of even higher Lorentz factors 
(F > 50, lAharonian et al]|2007l ). Assuming, that the mag- 
netic acceleration is efficient in these sources, and hence the 
observed F is close to ao, we obtain a characteristic value of 
ao ~ 10, which is much lower compared to GRBs. In prin- 
ciple, the value of ao can be higher near the black hole and 
then decrease downstream, e.g. as a result of some mass- 
loading process. However, at present we have no concrete 
evidence for this. Finally, the obser ved half-opening angle s 
of blazar jets are smaller, ~ l°-3° l|Pushkarev et al.ll2009l '). 

The corresponding rescaling of the results for GRB jets 
yields the light cylinder radius 

nc = 6 X 10- cm , (63) 

the distance to the fast magnetosonic surface 

'»"(S)"(lo^)"' 

and the shortest variability time scale required for establish- 
ing a steady-state super-fast-magnetosonic flow 

Since the timescale of strong central engine variability is 
unlikely to be shorter than the viscous timescale of inner 
accretion disk, 

tvis ~ 100i?,/c^ ~ 1.0 (j^^^ days , 

we conclude that, just like in the case of GRBs, the initial 
acceleration of AGN jets up to super-magnetosonic speeds 
is provided in a steady-stat e fashion. Moreover, th e recent 
observations of AGN jets l|Pushkarev et all |2009| ) clearly 
indicate that they satisfy the FOj < 1 condition of effec- 
tive st eady-state coUimation acceleration (|Komissarov et al] 
l2009al l. The observed decrease of half-opening angle with 
distance in M87 jet also supports the theory of coUimation 
acceleration l|Biretta et al.ll2002l : lGraciall2005l ). 

The distance at which one half of the electromagnetic 
energy is converted in the energy of bulk motion is now 

(see also iKomissarov et al.ll200'ii l. This scale is unresolved 
with modern VLBI systems. The recent numerical simula- 
tions show that the coUimation acceleration may continue a 
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bit beyond this point, reducing the magnetization down to 
cr ~ 0.4 within another decade of distance from the source 
(jKomissarov et al]|2007l ). This is still a relatively high mag- 
netization leading to relatively low efficiency of MHD shock 
dissipation. Additional impulsive accelerative can improve 
this. This time, however, when the impulsive mechanism 
switches on this is already the coasting regime. Indeed, 
the fast magnetosonic speed corresponding to a = 0.4 is 
/3ms ~ 0.5. Then instead of Ea. (l38|l the speed of the head 
rarefaction is given by 



2r2 + 1 



(67) 



and the length of the section of the shell affected by the 
rarefaction grows with time at the rate 



A/3h = /3 - A, 



2Vi 



+ r2 



The scale of transition to impulsive regime is now 



Rc 



For the typical parameters of AGN jets this reads 



Rc 



\wdj Vio 



pc. 



(68) 



(69) 



(70) 



Basically, since /3nis is mildly relativistic we have -Rcr.h ~ 
RcT,t, and since F ~ cto there, the two distances are also of 
the order of Rc ~ a-gctv. Thus, the theory predicts effective 
dissipation at internal shocks on the scales of ~ 1-10 pc, 
exactly the region where VLBI observations reveal bright 
super-luminal knots of AGN jets. This is the AGN counter- 
part of the prompt emission region of GRBs. 



6 SUMMARY AND CONCLUSIONS 

In this paper we investigated the properties of magnetic ac- 
celeration of relativistic impulsive ffows. As a first step, we 
focused on a relatively simple test case where a uniform 
cold and highly magnetized (ctq 2> 1) shell of initial width 
lo, whose back end leans against a conducting "wall" and 
whose head faces vacuum. The evolution of the ffow that 
develops in this test case splits into three distinct phases. 

The first phase can be described as a formation of a 
plasma pulse (or a moving shell). During this phase, which 
lasts for the time ~ to = lo/cms.o ~ lo/c, a self-similar rar- 
efaction wave develops at the interface with vacuum and 
travels towards the wall. At the end of this phase, the mean 
Lorentz factor of the outfiow is only (F) ~ and, apart 
from the very thin layer at the vacuum interface, the shell 
of plasma is still highly magnetized, with a mean magneti- 
zation parameter of (cr) ~ o-f/'^. 

The first phase ends when the rarefaction wave reaches 
the wall. At this point a secondary rarefaction wave forms 
that propagates from the wall into the back of the shell and 
decelerates the material that passes through it so that the 
shell quickly separates from the wall and moves away from it. 
During this second phase, the center of momentum Lorentz 
factor of the shell remains fairly constant (Fcm ~ ""o''^)- 
However, the leading part of the plasma shell, ahead of the 
secondary rarefaction, continues to accelerate at the same 



rate as in the self-similar solution. It contains most of the 
shell energy and its mean Lorentz factor grows as (F) oc t^''^ . 

At the end of the second phase, which lasts up to 
~ = o-gto, the magnetization of the shell drops down to 
cr ~ 1, one half of the electromagnetic energy is converted 
into the bulk motion kinetic energy of the plasma, and the 
growth of the mean Lorentz factor begins to saturate at 
(F) ~ (70 . Thus, the fiow enters a phase of coasting. During 
the coasting phase the pulse width grows faster, approach- 
ing I oc t. The decrease of the magnetization parameter also 
accelerates, approaching a oc t~^, and the pulse soon be- 
comes kinetic-energy dominated. This property of impulsive 
magnetic acceleration is most valuable in astrophysical con- 
text as the efficiency of relativistic MHD shock dissipation 
decreases dramatically with magnetization. In contrast to 
an impulsive flow, a steady-state magnetized jet either re- 
mains highly magnetized (u ^ 1) all the way, or approaches 
a ~ 1, depending on the efficiency of external colhmation. 
This implies at best only modest shock dissipation efficiency. 

In some cases of truly explosive phenomena, such as 
magnetar flares, our impulsive magnetic acceleration mech- 
anism can be solely responsible for the flow acceleration. 
In most other cases, such as GRB and AGN jets, strong 
variability of their central engines is not expected on time 
scales below the viscous time-scale of the inner accretion 
disc around a black hole, which powers relativistic outflow. 
This gives plenty of time to establish a quasi-steady super- 
fast-magnetosonic flow near the source where it is acceler- 
ated via the coUimation mechanism. The observed strong 
colhmation of these jets supports our conclusion that the 
colhmation mechanism plays a part in their acceleration. 
The impulsive acceleration mechanism comes in force fur- 
ther out, where an individual ejecta element starts being 
accelerated after the head rarefaction crosses it and creates 
conditions similar to those of our test case flow in phases 
two and three. The mean Lorentz factor of the shell, how- 
ever, starts increasing significantly above the value achieved 
by the quasi-steady coUimation acceleration only when the 
tail rarefaction wave crosses about half of the shell. Provided 
the central engine variability is sufficiently strong, so that 
the flow can be described as individual ejecta shells sepa- 
rated by long gaps, the impulsive acceleration mechanism 
can complete the acceleration process and produce kinetic 
energy dominated relativistic flows on astrophysically rele- 
vant distances from the central engine. For short GRBs this 
may still work well even if the ejecta effectively form a single 
uniform shell. 

Our analysis of GRBs show that a combination of the 
colhmation and impulsive mechanisms can accelerate GRB 
jets up to F ^ 10"^, as has been inferred recently for several 
bright GRBs detected by the Fermi Large Are a Telescope, 
for b oth long ( Abdo eta l. 2009a, b) and short (|Abdo et all 
l2010h duration GRBsEl Moreover, their jets can become 
kinetic energy dominated before the interaction with the 

We do note, however, that these lower limits on F from 
pair opacity are somewhat model dependent and a fully self 
consistent calculation appropriate for an internal shock ori- 
gin of the ga mma-ray emission gives limits that are a factor 
of ^ 3 lower l lGranot. Cohen- Tanugi fc do Couto e Silval |200S| : 
lAckermann et al.ll20icl) , F > lO^-S, which are significantly easier 
to satisfy. 
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interstellar or stellar wind gas begin to decelerate the 
ejecta at i?dec ~ lO^'^ — 10^^ cm. The dissipation at in- 
ternal shocks can become efficient on scales Rc ~ 
10^"^(cro/300)^(fi,/4ms) cm. The large variation of Lorentz 
factor at the coasting phase, AF ~ F, insures that the in- 
ternal shock will be strong and can dissipate and radiate of 
the order of ~ 10% or so of the flow kinetic energy, leading 
to a possibility of strong prompt emission. 

The AGN jets are likely to be accelerated up to their 
observed Lorentz factors already during the coUimation ac- 
celeration phase. However, the impulsive acceleration phase 
remains important, providing effective conversion of remain- 
ing electromagnetic energy and producing kinetic energy 
dominated flows. Our estimates show that efficient shock 
dissipation region, analogous to the prompt emission region 
of GRBs, is located around ~ 1 — 10 pc, where VLBl obser- 
vations reveal the presence of super- luminal "blobs" . 
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APPENDIX A: SELF-SIMILAR RAREFACTION WAVE IN PLANAR SYMMETRY 

The equations of relativistic MHD can be written as 

a^r^" = , a^F*^" = o , d^,{pu'') = o , 

fsee lKomisaroviri999l . and references therein), where 



T 



{phg + b )u'^u' + \Pg + 



9 



(Al) 



(A2) 



is the energy-momentum tensor. Here p, Wg — phg, pg, and u'^ — (u°,u) = r(l,'iT) are the fluid proper rest mass density, 
enthalpy density, pressure, and four- velocity, where F = (1 its Lorentz factor, g'^" is the metric tensor, and we 

use units where c — 1 for convenience. Furthermore, 



is the dual tensor of the electromagnetic field, and b'^ — {b'^, b) where 

r B + b°u B 



B-u = VB-v 



+ V{v- B)v , 



F 

is the four-vector of the magnetic field, which is defined as 

J, 1 P r^li 



(A3) 



(A4) 



(A5) 



where F'^^ is the electromagnetic tensor and rjap-yS is the Levi-Civita alternating tensor. In the fluid rest frame 6'' = (0, B) 
where B is the usual three- vector magnetic field, divided by -s/^tt, so that Wm — 2pm =6^. In general, B is measured in the 
lab frame. The three-vectors of the magnetic and electric fields in an arbitrary frame are given by 



B = F"' 



bu 



ub' 



E — b X u . 



Similar to classical MHD, the electric current is given by the second Maxwell equation, 

r = duF^" , 



(A6) 

(A7) 
0, Uub" = 0, 



where it also includes the displacement current (time derivatives of the electric field). Finally, \/ ■ B — diF*^' 
Uf.u'^ = -1. 

The RMHD equations simplify considerably under the assumption of a flat space-time, g^" = jy*"^ — diag(— 1, 1, 1, 1), 
and planar s ymmetry, i.e. that all quantities depend only on x and t in a Cartesian coordinate system (see, e.g., 
iGiacomazzo fc Rezzollal[2006h . 



d_ 

dt 



phV' -p- pV 

phr^v'' - b'^b=' 

phV'^v'> - b^b^ 

phr^v'' - b'^b'' 

B^ 



d_ 

dx 



p/iF^i;" - b^'b"" ~ pFi;" ' 
phT^v'^v'' +p~ b^'b'' 
phT'^v^v'' - 6"b^ 
phT'^v^v'' - b"6" 



(A8) 



Here we consider the even simpler case where = =Q and B^ — B^ = so that v = vx and B — By, i.e. u'^ — F(l, v, 0, 0) 
and b'^ = (0, 0, B/F, 0). Under these conditions, the RMHD equations further simplify to 



d_ 

dt 



( 



p/iF^ - P - pF 
B 



+ 



pVv 
phV'^v - pVv 
phVv' 
Bv 



-P 



= , 



(A9) 



and the the magnetic field in the fluid rest frame is given by B' = B/F, so that the equations for the evolution of pF and 
B = B'T are the same and B/pT — B'/p — const. Thus, we are left with three equations for three variables (e.g., p, F, and 
Pg), where we also need to assume an equation of state. In our notation 



hg 



P 



ha 



El 
p 



1 + 



7 Pa 
7- 1 p 



(AlO) 



where e = eint/p and eint = ep — ' 



Pg — p is the proper internal energy density of the fluid, while 7 is the adiabatic index 



of the fluid. The total pressure is given hy p — Pg + Pm = Pg + b^ /2. 

We are looking for rarefaction wave solutions, which are self similar, i.e. all quantities depend of x and t only through 
their ratio, which is deflned as the self-similar variable: £^ = x/t. In rarefaction waves the speciflc entropy, s, of every fluid 
element is conserved, and therefore = ds/dt — ds/dt + vds/dx. Since d/dx — (l/t)d/d^ and d/dt — —[(^/t)d/d^, this implies 
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{v ~ (,)s' — where a prime denotes a derivative with respect to ^ (s' = ds/d^), and therefore s' — and s = const (in general 
V ^ Therefore, the flow is isentropic, and we may simply relate the pressure to its value ahead of the rarefaction wave, 



p = poo-oP , (All) 



Po 



1 



Pa = Ps.o ( — ) = PoaoP^ , (A12) 



= "^Th'-'i^) ='»(''+7^""'"' + ''"''') • 

where p = p/po, tio = Vgfi/po, and ctq = Bq/pq is the magnetization parameter of the fluid ahead of the rarefaction wave 
(which is assumed to be at rest in the lab frame: To = 1). 

equation IA9I can be expressed in terms of the self-similar variable f as 

= {v-0{p' + P'i^^vv') + pv' , (A15) 
= {v-€,){phV'^)' + £,p + phV^v' , (A16) 

= (1 - Op' + {v- ^)phr^v' . (A17) 

Let Cs, ca, and Cms denote the sound speed, the Alfven speed, and the fast magnetosonic speed, respectively. We have 
— {\/hg){dpg/dp)s and — / ph, so that hg/h — 1 — c\ and 

cL = i(^|^) -cl + cUl-cl), (A18) 

which implies p' = c^^^^hp' . Therefore, Eq. HA17|I can be rewritten as 

= (1 - v^cLsfy + {v~ Op^^y' ■ (A19) 
equation. IA15I and lA19l implv c^s ~ [(« — 0/(1 ~ '^C)]^ and therefore 

c 

1 — ^4 



(A20) 



where the plus and minus signs correspond to rarefaction waves propagating to the left and right, respectively. This also 
implies 

1 =F t'Cms 

The velocities of the tail (where v — 0) and of the head (where Cms = 0) of the rarefaction wave are given by 

Ct = TCms , ?h = iWmax , (A22) 

where Umax ~ max is obtained at the head of the rarefaction wave. As expected, the tail of the rarefaction wave propagates 
into the fluid at rest at the fast magnetosonic speed, 
equation. IA19I and IA20I imply 

r^dv ± ^dp = =^ J±^-\n( ii^l ± f Hnill^dp' = const . (A23) 
p 2 \l-vj Jg p' 

Under our assumptions. 



^ / 'yapp ' ^ + ^op /.o/i^ 

Cms(p) = - / -r—^ --,-1 I ' (^24 

Y 1+ T^aop^ ^ +croP 

so that the integral in Eq. (|A23P can be calculated analytically in the simple cases where (Jq = {Bg = 0), or ao = {Pg,o = 0). 
In the first limit (tjo = 0, i.e. no magnetic field). 



2 



ilnfi±i^U_i_lnf^^±£i^ =const, (A25) 
l|Marti fc Mulleij|l994 ) so that 

/ , . ±— 2— 

' ' ' — const , (A26) 



1 — V J \ ^7 — 1 — c 
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and in the limit oo ^ 1 and 7 = 4/3 this imphes Fmax ^ 1 which is approximately given by, 



4700 


(7-1)-^/^ 1 


[(7-1) J 





2*^-'af (7 = 4/3) 



4ao = 4p3,o/po (7 = 2). 
In the second limit (ao ~ 0) we find 

J± = ^ In Q ± 2ArcSinh (^\/a^^ ~ \ (y^^) ^ In (^\/ aop + ^ aop + = const , 

so that 

and in the limit ctq S> 1 we have 

Tmax ^ SfTo • 



const 



(A27) 

(A28) 
(A29) 
(A30) 



It can be seen that the purely magnetic case, oo = 0, is equivalent to the pure hydrodynamic case, ctq = 0, for 7 = 2 and 
ao — >■ cro/2. In the more general case. 



J± = ^ In (^^) ± Hp) = ±^(1) = const , 



where 



7aop'>'-i + aop 



so that 



and 



v = ± 



exp[27(l)] - 1 
exp[27(l)] + 1 



exp[27(p)] - 1 
"exp[2/(p)] + l 



Tmax — 



exp[J(l)] ^ exp[2J(l)] + 1 
1 + Wmax 2exp[7(l)] 

exp[27'(p)] + 1 
2exp[7(p)] 



where 

7~(p) = 7(l)-7(p) 



dp / 7aop^-i + aop 
P \J 1 + ^aop^-"^ +aop 



(A31) 



(A32) 



(A33) 



(A34) 



(A35) 



APPENDIX B: THE AVERAGE LORENTZ FACTOR 

The maximal Lorentz factor Fmax is only asymptotically reached at the very head of the rarefaction wave, and only a small 
amount of material which carries a small fraction of the total energy has F ~ Fmax- Therefore, it makes sense to calculate 
some average value of the Lorentz factor, which would reflect better the Lorentz factor of the material that carries most of 
the energy. A natural definition is the weighted average over the energy. 



(r)i 



jTdE _J dx T°° F _ Jl' d^ T°° F 



JdE JdxTOO 



where 



rpOO 

1 = Po 



P + 



7-1 



aof + aof' 



~~f ao _2 
- aop^ - —P 



Another possible definition is the weighted average over the rest mass, 

/ VdM _ J dxT^p 



JdM JdxTp 



Itd^rp 



(Bl) 



(B2) 



(B3) 



We note that for a cold magnetized shell, at late times when almost all of the energy is in kinetic form and the magnetic 

energy becomes negligible, the enumerator approaches E/c^, and since the denominator is simply the rest mass M, then (F)m 
approaches E/Mc^ = 1 + ao/2. Alternative options to define a "typical" Lorentz factor are its average over space 
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Figure Bl. Right panel: four different estimates for the "typical" Lorenz factor, (F), within a rarefaction wave for the pure hydrody- 
namic case (co = 0) and for an adiabatic index of 7 = 4/3, as a function of ao = po/poc^- Various symbols show different weightings of 
r. Left panel: the same four estimates of the "typical" Lorenz factor within the rarefaction wave for the pure magnetic case (ao = 0), 

1/3 

as a function of ctq. The red solid line corresponding to (F) = CTq' has been added for reference. 



(B4) 



or its value at the point where there are equal energies on either side within the rarefaction wave in the lab frame (i.e. the 
"energy median" value), 



2:1/2 



dxT^ 



= \ r(Ci/2 



1/2 



€1/2 



(B5) 



1/2 



where Xmin = tit and a;max = tin (see Eq. IA22|) . 

Figure lBll shows these three estimates for the typical Lorentz factor within the rarefaction wave, for the pure hydrodynamic 
case ((To = 0; le]t panel) and for the pure magnetic case (ao — 0; right panel). In the pure hydrodynamics case the typical 
Lorentz factor of the material in the rarefaction wave is only mildly relativistic even in the limit of ao 3> 1, where it approaches 
a constant value, while Fmax rapidly increases with ao (see Eg. IA27] ). In the purely magnetic case, we find that the typical 



value of the Lorentz factor within the rarefaction wave is (F) 

wave is Fmax ~ 2(70. 



1/3 



while its maximal value at the head of the rarefaction 



APPENDIX C: ANALYTIC DERIVATIONS FOR THE RAREFACTION WAVES 
CI Explicit solution for the original self-similar rarefaction wave 
Equation ((5| can be written as 



and Eq. ((6]) as 



where Sx is defined via 



1 + X 



This allows us to find all flow variables as explicit functions of ^: 



X : 



u = Tv = 



(^c„,o(?,)"/^ - ] 



5c„., = 



^2/3 



1 + (T 



^4/3 ^-2/3 



54/3 ^-2/3 



+ 1 



(7_p_B'_ B _ 1 I S. 
CTo po Bo TBo 4(70 



2/3 



0^ 5c,„,,o , 



^2/3 



(CI) 

(C2) 



(C3) 
(C4) 

(C5) 
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(this result is due to lLvutikovll2010al . with a smaU correction in Eq. (|C5|I ) . The extent of the rarefaction wave is given by the 
conditions (5c„, = 1 for the right front and 5„ = 1 for the left front. They yield (5~^_. ^ < <5f < <5c„-,s o- 



C2 Motion of the head of the secondary rarefaction wave 

The overall impression created by Fig. [2] is that of an effective separation of the shell from the wall. The region between 
X = 18Zo and x — 20Zo contains most of the total mass and energy of the initial solution. This somewhat surprising result can 
be verified in a different way. In the fiuid frame the front of secondary rarefaction moves with the local magnetosonic speed. 
In the lab frame this corresponds to 



dx. 



cS/3 c2/3 



(C6) 



where ii(^) and Cms(C) describe the self-similar solution for the initial rarefaction (Eqs. [C3) . |C4| V Noting that d£,,/dt — 
iJdt = {d5lJdi,){di,/dt) = {d^,/dt){5l^ 



[{dx,/dt) - and dSjjdt = {dSlJd^,){d^,/dt) = (d^,/dt)(5? + if /2, equation can be rewritten as 



dSl 
dint 



sL + 1 



(C + 1)^ 



which has the solution (|Lvutikovll2010al ) 



1 - 



2/3 



-3/2 



where 5^^ < 



and 



For (To ^ 1 we find that to,* 



to,. (^cU-lf^^ 

to <5ct..o('5L=.o + l) 
Iq so that at t S> to we have (5? r 



4ag/^(l + ao)^/" 
2/(1 -e.) > 1 and 



t 



1 



1 - ^* . 
For (Jo 3> 1 we also have 5l 



-3/2 



1 + 



t 



16agto 



-2/3' 



3/2 



2to/t 



t < 16croto 



' 4(To, so that for to <^ 16(Toto, where 5'i 



Ufo/ 



r. = (i-/?i) 



-1/2 



I 1/8(72 

t/to, equations (|C3|) and (jlOl) imply 



to 



(C7) 



(C8) 



(C9) 



(CIO) 



(Cll) 



Note that F* is the Lorentz factor of the motion of the head of the secondary rarefaction wave, while F(^,) is the Lorentz 
factor of the fluid at that location. 



C3 Analytic calculation of physical quantities at ^ > ^*(t) 

It is instructive to calculate the values of relevant physical quantities in the region between the head of the secondary 
rarefaction wave and the vacuum interface, which corresponds to ^*(t) < £, < /3max. In particular, it can help verify that this 
region contains most of the energy and rest mass in the flow during the magnetic rocket acceleration phase (at t <^ tc = o-Qto). 
Using equations (|C3|) and (|C5I) . at a given time t we have dx = td^ so that the rest mass per unit area at ^ > 5, is given by 

2 



M[> Ut)] 



c(t) 



dxr{x)p{x) 



pot 

8(70 



d^ S, 



2/3 r2/3 




where 



/3n 



2^(7o(l +(7o) 



(C12) 



(C13) 



1 4- 2(70 

is the maximal fiuid velocity, which is obtained at the vacuum interface, while ^*(t) is given implicitly by equation (|C8I) . One 

x2 



can change variables of integration to 5^ and then to y 



2dSl 



S^^^ , using the relations 
Gy^dy 



(Sj + 1)2 (y3 + 1)2 ' 



(C14) 



and use the simple expression for the initial mass. Mo = polo = poto[(7o/(l + cfq)]^^^, to obtain the following expression for 
the fraction of the initial rest mass at ^ > £,*{t), 
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4(7, 



3/2 
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[1] 



(C15) 



where ymin(i) = '^j^l't) ~ ^^i^s o- This result can be written more explicitly and simplified using Eqs. HC8|) and (|C9|) . 



Af[>C*(t)]_ 4(1 + ao) 



Mo 



(Vl + cro + ^)' 



1 + 



c2 X 1/3 



(C,„ + i)(^) 



1 t < cToio , 

2(Joto/t t 2> (7oto , 



(C16) 



where the last asymptotic values are valid for cro 2> 1. 

Using Eq. (|C10|) . one can calculate the fractional change in the width of the region between the secondary rarefaction 
wave and the vacuum rarefaction, A,(t) — t[(,h — where — /Jmax, 



A, (to) toKh-Ct] 



1 + 



V 1 6(7^*0 



2/3' 



3/2 



1 + i3/2){t/16a^tof/^ t < IGaito 



Wa^to y (3/2)(t/16<7gfo)'/' t » 16agto 



(C17) 



where = — Cms.o ~ — [cro/(l + ao)]^^^, A, (to) ~ toCms.o iXtcro ~ ^'^'^ asymptotic values are valid for tJo ^ 1. 
The kinetic energy, Skin = / da::r(r — l)p, is given by 



gkin[> C*(^)l ^ 3^/TT^ / t 
Mo 



dy-, 



Vl + cro t 

8, 



ay 



«y + — Uy + 2 2 + ^ 



ay 



y a- 



2a + 6a^ + 2a^ - (4a^ + 6a^ + a**)y + (1 + 6a^ + 4a'^)t/^ 



f{y) = 



+ - 



1 - 2a2 - 2aS + a* 



a4(l + y3) 
^/3 



■ arctan 



, , , (a^ + l)3(a2 - 1) , 
+ 3 Hy) + In 



1 + y 



(CIS) 



(C19) 



The electromagnetic energy, -Bem = / da;po<7o(r'^ — i)(p//9o)^, can be calculated in a similar way, 



EEM[>^4t)] _ 2 eem[> e*ft)] ^ syrr^ / 1 \ r' ^ 



-Eem.o 



CTQ Mo 



(y3 + 1)2 



2 2, 1 

a y + 



a^y^ I \ y 



o- ^ y 



5(y) = 4--- 



16cr,^ 

y a^ 1 + Sa*' + a^^ - Ta-'y - 4a^°y + Aa^y'^ + 7a**y^ 



a y 
4(1 - 2a2 - 2a'' + a*) 



arctan 



3a6(l +y3) 

2y-l 



, , 4(a2 + l)^(a2 - 1) , 
41n(y) - ^ '- In 



1 + y 



^l + y2 -y 



(C20) 



(C21) 



The total energy (including rest energ y) is simply / T^''°dx = M + iJkin + -Eem- This can be used for the normalization when 
calculating the average values of quantities weighed by 



rpQO , N PO 



4(ay+-l) ^-2+-^| + (aV + 



1 \ a 



— -2+4. 



We find that 



E{>i,{t)\ 2E{>i,{t)\ VTT^ / t \ (a^ -y^i„)3 [y^i„(l + a6)_|_3„2(.^^^2^2^.^)j 



So (2 + ao)Mo 16(2 + cro)cr;^ 



a6y,^i„(y^i„ + 1) 



(C22) 



(C23) 



2/3 2/3 

and the same holds for the energy above some ^ > ^*(t) where ymin(t) = \t) replaced by y(t) = 5^ . In order to calculate 
the mean Lorentz F factor or magnetization a at ^ > ^, (denoted by (F)* and (cr)*, respectively) one needs to calculate the 
following integrals: 



1 

Mo 



dxr°°F: 



3/2 



(I;) [/r(«')-/r(ymin(t))] , 



(C24) 



My) 



3(a2-y4) 9a*(l + a«) + a2(l-25a«+ai2)y-(l-25a«+ai2)y2 3(a'' - 1) , ,^ , 3, 3(a« - 2) , , , 

■ + y ) + 64^3 ln(y) 



128ay2 



64a7(l + y3) 



64a3 
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64^3 c 



V3 



l + y 



192a'^ 



Mo 



a:vac(t) 



(C25) 
(C26) 



3{a^ - Aa'^y + Aa^y^ - y'^) 2a^(l - 250*^ + a^^) - (1 - 340*^ - 8a")y - a^(8 + 34a'^ - a^^)y^ 3 
•^"^^^ " 256^V ^ 128a8(l + y3) ^'^^S'^ 



1 + 10a* - lea"^ - W + lOa^^ + a^^ /2u-l\ (a^ + 1)^(1 - Sa^ + 5a* - a«) , 

■ arctan | ^— | + — — „„ , „ In 



128\/3a8 



%/3 



l + y 



384a« 



(C27) 



APPENDIX D: GENERALIZING TO A SPHERICAL FLOW 

Here we consider the case of cold (pg = 0) radial flow. We assume that the flow is one-dimensional and the magnetic fleld is 
perpendicular to the radial direction. Obviously, this is not fully self-consistent, but this is a reasonable approximation for an 
equatorial wedge. More accurate two dimensional treatments are saved for future works. 
The basic equations for a one dimensional flow in spherical symmetry are 



dt 



dt{pr) + —dr{rprv) 

dtibV) + -drirbTv) 
r 

+ idr[r\p + b')r\] 



(continuity) , 
(magnetic field) 

(energy) , 
(momentum) , 



where v is the velocity, b = B/\/4TTr, and B is the (azimuthal) magnetic field as measured in the source frame. 
One can introduce new variables, b, p and x as follows 



-2- I -17 

p — r p , = r 
which upon substitution into Eqs. (|DHD3[) yield 

dtipT) + d^ipvv) = 0, 

dtibV) + d:,(bTv) = 0, 

52 



dt 



{p+b')V' 



dt [{p + b')v\] 



+ d. [{p + b^)r^v] 
(p + -,^)rV + 5 



, 



= 



(Dl) 
(D2) 

(D3) 
(D4) 

(D5) 

(D6) 
(D7) 

(D8) 
(D9) 



These equations are identical to those of plane cold {pg — 0) fiow. Therefore, all the results obtained for the planar case, 
including the self-similar solution, can be utilized in the spherical case. 

After the substitution p — >■ P/2, Eqs. HD6|I . (|D8|) . and HD9P also become identical to those of unmagnetized plasma. From 
Eqs. (|D6P and (|D7P it follows that b/p = const, and p ex. p^ . Thus, the ratio of specific heats for this plasma is 7 = 2. 



APPENDIX E: SELF-SIMILAR SOLUTION FOR THE SHELL'S TAIL. 



Our starting point is Eqs. (|D6|) - (|D9p . which are valid for both the planar and spherical (after the substitution (|D5p ) cases. 
Equations (|D6|) and (|D7|) imply 

d_fl 
dt\p 



, (El) 



that is b/p remains constant for each fiuid element, and is determined by the initial conditions. 

Just like in Appendix]^ we introduce the self-similar variable ^ = x/t, but this time we seek solutions of the form 

v = V{0, P^t-FiO, b^t-GiO (E2) 



(Since b/p remains constant for each fluid element, b and p must have the same temporal scaling.). Using Eqs. (|D6|) and (|D7|) . 
Eq. (|D9|) can be reduced to 
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/ dvv drv\ 2 dv idb^ 



Substituting the expressions HE2p into this equation we obtain 

t^FTiV - C)(rV)' + t^^iGG' - G^T^iV') = 0. (E4) 

This equation is satisfied for any t only in the following two cases. First, if q = - this is the case of simple rarefaction wave 
analyzed in Appendix [A] Second, if 

V^i, (E5) 
and 

GG' - G'^V'^ey' = . (E6) 
Integrating the last equation (after substitution of equation |E5) into it) we find that 

G(C) = = AV{i) . (E7) 

Substitution of expressions (|E2I) into Eqs. HD6|I and HD7|I leads to 

(l + a)F = 0, (l + a)G = 0, (E8) 

which are satisfied when a = —1. Function -F(^), however, remains undefined. 

If the initial solution is uniform, as in the planar case considered in § [2l with po and bo being the initial rest mass density 
and magnetic field respectively, then from Eq. (jEl[) we obtain 

p=f&, (E9) 
bo 

and thus, 

nC) = fG(e) = f^£=. (ElO) 

bo bo _ ^2 

Thus, a self-similar solution of the required form does exist. In planar geometry, this is 
and in the spherical geometry 

In both cases the magnetization parameter of fiuid elements decreases linearly with time: 
J? 

a^ — ^ — xt-\ (E13) 
P P 

Clearly, this solution is only applicable for ^ < 1. Moreover, it cannot be simply truncated at some large ^ and continued 
with vacuum. Instead, it should terminate at a shock or smoothly transform into a non-self-similar fiow. For our simple test 
problem, this solution cannot become asymptotically valid up to f = /3max given by Eq. HC13|) at very late times {t 2> tc), 
since it cannot simultaneously satisfy the global conservation of energy and rest mass for an initially uniform shell. This 
implies a limited region of applicability, in the tail of the fiow but not at its head0 General considerations, however, show 
that some similar scalings still apply at the head of the flow at late times: a oc while p oc in planar (spherical) 

geometry. Moreover, v = ^ at asymptotic late times when the magnetic pressure becomes dynamically unimportant and each 
fiuid element moves ballistically. 

Finally, we calculate the scalings of the total kinetic -Ekin(i, Ci) Cz) and magnetic -B_g(t, Cii Ca) energy of a fiuid element 
bounded by xi = ^it and X2 — (,2t. Both in planar geometry, where dV oc dx = td^, and in spherical geometry, where 
dV oc r^dr — t'^^^d^, we obtain 



£kin(t,Ci,6) = j dVpTiV ^l)oit\ (E14) 
Si3(t,Ci,6) = /'dl^&'(r"-i) oci-i (E15) 



Likely, once cr{£^*) drops below unity the flow can no longer efficiently rearrange itself from the original self-similar structure (described 
in Appendix [A]| to the new one (described here), as the secondary rarefaction wave becomes weak). 
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APPENDIX F: ACCELERATION AFTER THE SEPARATION FROM THE WALL 
Fl Mechanical analogy: two masses and a spring 

At t < to the acceleration of the shell as a whole occurs mainly because the back end of the shell is pushing against the wall, 
and therefore this mode of acceleration can remain effective only as long as the shell remains in causal contact with the wall. 
Therefore, the initial shell crossing by the rarefaction wave accelerates the shell up to (F) ~ , and soon after the shell 
can no longer effectively push against the wall, as the magnetic pressure at the wall drops dramatically, and the subsequent 
change in the total momentum P of the shell due to the force F exerted on it by the wall {dP — Fdt) becomes negligible. 

It is a somewhat surprising result that after the shell separates from the wall its mean Lorentz factor continues to increase 
with time despite the apparent lack of any external force. This can be understood as follows. The total energy and momentum 
of the shell are indeed conserved in the lack of an external force (or energy losses or gains). However, the shell expands 
under its own pressure, and develops a considerable relative velocity between its leading and trailing edges. In its center of 
mass frame the energy and momentum of the front and back ends of the shell are comparable. However, if the expansion is 
relativistic in the center of mass (or comoving) frame then in the lab frame the energy and momentum of the leading part 
are much larger than those of the trailing part, and thus the leading part dominates the total energy and the Lorentz factor 
when averaged over the energy in the lab frame. 

This may be illustrated by the following simple example. Consider two identical masses m moving together with a 
compressed ideal massless spring between them, with potential energy i?pot in its own rest frame (S, , which is also the 
rest frame of the two masses, hereafter the comoving frame). The energy of the system in its own (comoving) rest frame 
is E' — 2m(? + E-pot, and in a frame where this system is moving at a Lorentz factor F = (1 — Z?^)^^/^ in the positive 
x-direction (hereafter, the lab frame), its energy is E — TE' = V{2m(? + -Epot) and its momentum is P^ = VjiE' /c = 
F/3(2mc^ + Epot)/c = /3E/c in the x-direction (while Py — P^ = 0). Then the spring is released and all of its potential energy 
is converted to kinetic energy of the two masses, which in the comoving frame now move at a Lorentz factor F* = (1 — /3^)~^/^ 
such that -Bpot ~ 2(F, — l)m(? and F* — E' /2mc? , in the positive and negative x'-directions, respectively (the two masses 
are thus denoted by subscripts '+' and ' — ' accordingly). In the comoving frame their energy-momentum 4-vectors read 
y'^M _ r,(l, ±/3,, 0, 0), and a simple Lorentz transformation shows that in the lab frame 

< = [FF*(l±/3/3.),FF.(/3±/3,),0,0] , (Fl) 

which indeed satisfies E ^ E+ + E- ^ mc^{u% + u'L) = 2FF, mc^ = TE' and P:^ = mc{u\ +u\) = fiEjc (and Py ^ P^, ^ 0), 
as it should, while F± — E±/m(? = FF*(1 ± /3/3,). Thus, the ratios of the energy and momentum of the two masses, and the 
fractions of the total energy and momentum that each mass holds are given by 

g+ _ F+ _ 1 4- PP, P^+ _ V+P+ _ P + P, E± _ 1± PP, P^± _ P± P, 

E^ F_ 1-/3/3. ' P^- r_/3_ P~P, ' E 2 ' 2/3 

For F, F, » 1 we have 

E- 1/1 1 \ , F4 



(F2) 



E+ 



+ ^-2F.»l. (F3) 



Therefore, almost all of the energy in the lab frame is in the leading mass (or leading part of the shell), which has greatly 
increased its Lorentz factor. For F, = F the only energy left in the trailing mass (or trailing part of the shell) is its rest mass 
energy and all of the potential energy is converted into the kinetic energy of the leading mass, which in this case also carries 
all of the momentum. Thus we can see that the leading mass, which constitutes one half of original rest mass, ends up with 
almost all of the energy and with a much higher Lorentz factor than what it started with. Going back to our magnetized shell, 
the potential energy in a "spring" is the analog of the magnetic energy in the shell, and similarly to the mechanical analog 
eventually most of the energy ends up in a good fraction of the original rest mass, that can reach a very high Lorentz factor 
(much larger than the initial Lorentz factor of the shell). 

F2 Evolution of (F) after the separation: an alternative derivation of the scaling (F) oc t^^'^ 

Here we follow the mean shell parameters but drop '()' in the notation, for simplicity. Let us consider a planaiF^ shell initially 
(at lab frame time t — 0) at rest in some rest frame So, which we refer to as the lab frame, in which it has a width Zq, magnetic 
field -Bo, rest-mass density po, magnetization ctq = Bq/4,-kpqc? ^ 1, energy Eq, and no (or negligible) thermal pressure (for 
simplicity). This shell can either be leaning against a wall to one end (at x = — /o), or be half of an unbounded shell (initially 
occupying — 2Zo < x <Q). We have shown that initially the shell expands due to the passage of a self-similar rarefaction wave, 
which crosses the shell over a time to ~ lo/c, and is accelerated to a typical Lorentz factor of Fi ~ (^^^ ■ 

Now, even though the shell is no longer perfectly uniform as in our initial configuration, we consider the part of the shell 
that carries most of its energy (as measured in 5*0), which is expected to be roughly uniform (the relevant physical quantities 

not changing by more than factors of order unity within that region), and make the analogy between it in its own rest frame, 

1/3 

5*1 (which moves at a Lorentz factor ~ Fi,o ~ Oq relative to So), and our original configuration (that was quantified in So)- 
The same same reasoning essentially also holds for a spherical geometry, as demonstrated in Appendix [P] 
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Even though this analogy is not perfect, we still expect a similar qualitative behavior, and a similar quantitative behavior up 
to factors of order unity (which we discard here, as we are interested only in the relevant scaling laws). 

One difference, however, is that for reasonably smooth initial conditions we no longer have a strong rarefaction wave 
crossing the shell, which eventually splits it in two (as for a perfectly uniform shell with sharp edges surrounded by vacuum 
on both sides, where the two rarefaction waves from both sides meet and are secondary - our "wall" for a one-sided shell). 
Thus, the shell is basically the smooth peak of the lab frame energy density (which scales as when cr S> 1; see the profile 
of B at t = 70to in Fig. [5]). The shell still significantly expands in its own rest frame, reaching speeds of order of its fast 
magnetosonic speed on its fast magnetosonic (or light) crossing time. However, this spreading is smooth and continuous, and 
the shell does not split in two, but instead it remains a smooth peak of the lab frame energy density. Material ahead of the 
peak in the pressure (at the front of the shell) is accelerated, while material behind this peak (at the back of the shell) is 
decelerated (by the pressure gradient, in both cases). 

Nonetheless, it is instructive to divide this process into discrete steps or phases. The approximate initial conditions of the 
"second phase" in the evolution of the shell, as expressed in frame Si (using a subscript "1" for all the relevant quantities, 
when measured in this frame) are 

Bo po 2/3 Eo 
— ' Tj3 ' <7i ~ o-Q , El ^ . 



In frame 5*1 , the two sides of the shell are expected to accelerate in opposite directions, and develop velocities of the order of 



the shell's magnetosonic speed, r2,i ~ o" 



1/3 



2/9 
^0 . 



on the shell's magnetosonic (or light) crossing time, ti ~ Zi/c ~ CQ'^^to 



(as measured in frame Si) 
.1/3, 



The same scalings as before should approximately hold here as well. 



Bi 



h ~ CTj ii ~ (Jo lo, B2 jyg gyg , P2 J73 ^ , ~ cr^ ~ , E2 

We note that in frame So (i.e. in the lab frame) almost all of the energy ends up in the front part of the shell (which was 
accelerated in the direction of motion of Si relative to So), whose rest frame S2 moves at a Lorentz factor r2,o ~ ri,or2.i ~ (Jq^^ 
relative to So, and has an energy of £2,0 ~ ^2.oE2 ~ Eo as measured in So. The back part of the shell has a Lorentz factor 



El 



Eo 



1/3 



5/9 



of ~ ri,o/r2,i ~ (Tq in frame So, while its energy in frame So is ~ E2 ~ Eo/(Jq , i.e. a factor of ~ CTq ^ 1 smaller 

than that of the forward half of the shell, so it can be safely discarded as we are interested in the part of the shell that carries 

2/3 

most of the energy in frame So- In frame So, the the second phase of acceleration takes a time ti,o ~ Ti.oti ~ ctq to- During 



this time F increased from Fi,o 



1/3 + r 
• OrJ to I 2,0 



5/9 



r2,o 



2/9 



i.e. by a factor of ■ 

1/3 



2/9 



implying 



F oc t 



1/3 



(F4) 



Similarly, recursively repeating the same procedure n times, it can be shown that 



1/3 



l(7„ 



l-{2/3)" 



t/o. 



where C7„ = (Z„, i„, F„,o, -B„ ^, p„ ^, cr„ ^, _E„ ^) and Fo,o = 1. Noting that t„,o ~ F„,otn ~ i^/io, this implies that 



log g-n 

log GO 



and 



(I)' 



logF„,o 

log CTo 



l°g (^) 



log(Fn,o/Fi,o) 
log(f„_i,o/ta) 



log Go 

_ 1 

" 3 



1 - 



(I)' 



■°.^(fer) 

log fjQ 



n \ n — 1 



F oc t 



1/3 



(F5) 



It can also be seen that in the limil|__| of n ^ 1, F„,o — >■ o-q, i.e. the Lorentz factor approaches its asymptotic value that 
is achieved when cr ~ 1 at tc ~ ctq^o- This implies a transition radius to the coasting phase of ~ g'qIo ~ T'\lo, where 
Fc = F(i?c) ~ Go- Up until this time the shell width in the lab frame (So) remains approximately constant, Z„,o ~ 'n/F„,o ~ ^o- 
This derivation relies on the fact that the shell, which represents the leading part of the flow, carries most of the total 
energy and rest mass during the impulsive acceleration phase, so that its energy and rest mass are practically constant. As 
we have seen in §[22] this condition is satisfled for at least as long as the shell remains highly magnetized, g^\. 



APPENDIX G: THE CENTER OF MOMENTUM FRAME AND CALCULATIONS IN DIFFERENT 
FRAMES OF REFERENCE 

Let us consider a frame of reference Si moving at a dimensionless velocity /3i in the positive x-direction relative to the lab 
frame , and let us denote quantities measured in this frame with a prime. Also, let {ta, Xa) = (ta, a^a) = (0, 0) correspond 



In practice the approximation that cr 1 breaks down at it = CTj, 



' a few, after ^ Iog(Iog CTniin/ log ""o)/ log(2/3) steps, or ^ 4 — 6 



for (TO ~ 10 and (7^ 



1.8 — 3.9, so even for cro ^ 1 the shell becomes kinetic-energy dominated within a rather small number of steps. 



© 0000 RAS, MNRAS 000, 000-000 



30 Granot et al. 



the event 'a' of exposing the original front of the magnetized shell at rest to the vacuum, i.e. the onset of motion of the shell 
material. It is easy to show that the event '6' of the original rarefaction wave hitting the wall corresponds to {tb, Xb) = {to, —lo), 
where to = lo/cms,o = ?o[(l + cro)/cro]^^^, and {t'l,, x'^) = tori(l + /3iCms,o, —fii — Cms.o,). Initially (at t' = 0) the sheU width, 
density and rest-mass are I'q — lo/Ti, pQ — Fipo and AIq = polo = p'oi'o, respectively, so that 



T'' 



-po(i + (Jo)riui 



j/ rjil fix 

Lol 



-Mo(l + (Jo)ui 



(Gl) 



Now, between t' = Q and t' = t'^ = T -\- /3iCms,o)^o[(l + o-o)/o-oY^'^ the momentum increases due to the external force exerted 
by the wall, dP' /dt' = poo"o/2 so that at t' — t'l, when the rarefaction wave reaches the wall and the wall is removed (so that 
from that point on there are no external forces and the total momentum P' and energy E' remain constant) we have 



p'{t' >ti) = H + t'b 



-A/o(l + o-o)Mi+Moriy 



1 + cro 



(G2) 



We are interested in the center of momentum frame in which by definition P' — 0. Note that here, in contrast to the previous 
subsection, we evaluate the total momentum simultaneously in this frame, rather than in the lab frame. According to Eq. (|G2|) . 
this corresponds to 



/3cM = 



2 + (TO 



FcM = 



2 + cro 
V4 + 3(70 



(G3) 



At the CM frame (which moves at j3i = /3cm relative to the lab frame), P'{t' > t'^) = 0, so that the total momentum remains 
zero from the time when the shell separates from the wall. In this sense, the shell as a whole simply does not accelerate in 
the CM frame at t' > t'^. Moreover, the energies of the front part and the back part are comparable in this frame. The total 
energy at t' = is 



n/.OO 



Po 



r?(l + (7o) 



0£ 

2 



E'o = I'oT 



I rpl. 00 



Mo 



ri(l + ao) 



(To 

2ri 



(G4) 



and between t' = and t' — t'^ it decreases due to the negative work performed on it by the receding wall at its back, 
dE' /dt' ^ F ■ v = -/3ipo(To/2, so that 



E\t'>t'b) = E'o-t'bP,f^^M, 



-r 0) 



-Mo — ui 



l + CTQ 
(TO 



= AfoFi 



1 + 



(To ~ (Toj 



.(G5) 



For /3i = /3i_cM given by Eq. (|G3|) this reduces to 



E'cuit' > t'b) = ^^i±^Mo = 



Fl.CM 



(G6) 



The self-similar solution describing the original rarefaction wave can be expressed in a rest frame S" moving at a velocity 
13-w in the negative x-direction relative to the lab frame (where the wall is at rest), so that in this frame the wall is moving at 
a speed /3m in the positive a;'-direction. This implies 



* " t' " 1 + /3»C 



r r2/3 cl/3 r2/3 



^ t 1 - /3»e 



v 



Ov — 



2(5, 



2/3 rl/3r2/3 



.4/3 .2/3.4/3 
^2/3 51/3^2/3 



(G7) 
(G8) 



o 

(To 




(G9) 



(GIO) 



Using this result one can rewrite the integrals in § IC3l in terms of quantities in frame S' . In particular, they retain the same 
form up to the following simple substitutions: 



t^t 



y 



I .2/3 

y =5J, 



I .2/3 .1/3 ,,\ ^ I /,'\ jc2/3 

^ =''cL,o"™ ymin(t) ymm(i ) = 0^/(t/) 



(Gil) 



Finally, £,',{t') can either be computed either using the solution for (,*{t) and the relation Cl(f') ~ [C*(i) + /3™]/[l + l3w(,*{t)], 
or by directly generalizing the derivation from § IC21 as follows: 

^8/3 ^4/3x2/3 



dx: v'iO+CmsiC) 



dt' 



dSl 
d\nV 



= Si, + 1 - 



(G12) 
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Figure Gl. Left panel: The ratio of electromagnetic to kinetic energy in a frame moving at fim in the negative I'-direction relative to 
the lab frame (i.e. the rest frame of the initial shell and the wall) at the time ip when the rarefaction wave reaches the wall, as a function 

of 5^, for fTQ = 10^, 10^, 10*. Indicated by vertical lines are the lab frame and the CM frame, where this ratio is ~ "'q'^^ ^'^'^ ~ "^0^^' 
respectively (indicated by short horizontal lines). Middle panel: the same for the CM frame as a function of ctq, using Eq. IIGI8I I: this 
ratio approaches cFq^'"^ (dashed red line) for ctq 2> 1. Right panel: the evolution of the same ratio as a function of t' /t'^ in the CM frame. 



which has the solution 



4 



2/3' 



52 



-3/2 



3_0 < < 5u,5c„ 



s,0 



:,o and 
1)3/2 



to 



{St 



yi/2 



where we use the notation tg = fj, and the relation 

to 



to,, ^ 

SwtO Sw5t^^ g{Sc^^ g + 1) 

(sLs.o + si) 



Aal^\l + aoy^^ 



For (To ^ 1 (and 5iu <C (<5c„3 □) ^^d that t'o^, 

-3/2 



t'n 



1 



1 



t'o so that at t' 3> t'o we have 5|/ a 

2t'o/t' 



2/(1 > 1 and 

t' < Wa^Sit'o 



1 + 



'2/3 



t' > IGaiSlt'o 



(G13) 



(G14) 



(G15) 



(G16) 



so that t'c/tQ ~ (TqSI or fc/fo ~ ao for 5,1, = 5ii,,CM ~ o"q Since tp « to/S^ this imphes tj. ~ (To^^jfc 
It is useful to calculate the ratio of magnetic to kinetic energies, 



EkM[>e*it')] _ 3 

EiJ>a{t')] 4 
For the CM frame we have 



9[{a'f]^9[yLin{t')] 



-E'em(^o) 

KM 



2/3r-2/3>, 



2/3 r-2/3x 



/3u 



-/3c 



•\/o-o(l + (To) 



>,CM 



2 + (To - ^0-0(1 + 0-0) 3 



(2rc 



3(5: 



2 + ao + v/c^oll + fTo) 4(To 
We also generalize the result for the total energy in the flow, 

E'[>e*{t')]_ ^/TT^ ft'\ (a"-ymi„)' Kin(l + «") + 3a'2(l + a'Vnin)] 



E'it'o) 



16(2 



\ 3/2 

- f^ojf^o 



, , V go(i+°-o) 

J- + Piu 2+<T„ 



(G17) 



(G18) 



(G19) 



and the same holds for the energy above some ^' > where y'^init') ~ S^/^^i-^ is replaced by y'{t') = 5^!^. One can verify 

that this ratio is indeed 1 at = t'g when y^in = {Sw/Sc^^ 0)^^^ ~ S^ja! . 

Figure [Gil shows SEMZ-E'kin first at t!^ as a function of the velocity of the primed frame of reference {left panel), then at 
to for the CM frame as a function of ctq {middle panel), and finally for the CM frame as a function of time t'/to {right panel). 

Figure [G2l shows various estimates for the typical Lorentz factor in the lab frame as a function of time {left panel), along 
with the cumulative distribution of energy as a function of the flow Lorentz factor (or velocity) at different times, both in the 
lab frame {middle panel) and in the CM frame {right panel). The first two panels help to quantitatively address an important 
point that has been raised in §|4]in the discussion around Eqs. (|34|) and (|35|) . Taking the energy weighted average over F, 
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Figure G2. Left panel: various estimates for the Lorentz factor as a function of time for ctq = 10'^; V f denotes the Lorentz factor 
below which there is a fraction / of the total energy. For example, Fq.s = (r)^ j^ed is the value for which there is equal energy in faster 
and slower material. Also shown for reference are r(5*) and (r}« (which is (r)g calculated over the region g > Note that V j can 

be calculated analytically only as long as V j > r(^»), since otherwise we must include in the calculation the region ^ < which is 

not described by the self-similar solution. Middle panel: the corresponding cumulative distribution of the fraction of energy in the flow 
as a function of T for log]^o(t/to) = 0, 0.5, 1,..., 8 (red lines). The lower-right boundary (thin solid black curve) corresponds to r(^»). 
The horizontal blue line corresponds to Fq.s = (r)^ mcdi ^■nd the vertical blue line shows the value of F where it meets with the curve for 
F(5*): Fo.5 = F(^«) (TO- Right panel: the corresponding cumulative distribution calculated in the center of momentum (CM) frame 
at log^o(*'/*o) ~ 0-25, 0.5, 4.5 (red lines). The bottom branch for F'(^J,) (thick solid black line) and t' = t'^ corresponds to material 
with a negative velocity in the CM frame, which initially carries almost half of the total energy in this frame. 



(r)_B, is not a unique choice, and averaging over the 4-velocity u, {u)e, would give a similar result. However, as shown in §|3] 
using {I3)e would give a very different result, where at late times {P)e — ^ Pcm corresponding to Fcm = (1 — /?cm)~^^^ ~ '^o^'^ 
while {r)_E ~ {'^)e ~ f7"o- Fortunately, we can also estimate the typical value of F of the material that carries most of the 
energy in the lab frame without having to perform any averaging, thus avoiding the need to choose a specific function of the 
flow velocity to average over. The left panel of Fig. IG2I shows the median value of F, {F)_B,mGd = Fq.s {thick soltd blue line) 
according to Eq. (|B5p . as well as the values of F below which there is a fraction 0.2 (Fo.g; green line) or 0.8 (Fo.s; cyan line) 
of the total energy. The middle panel shows the corresponding cumulative distribution of the fraction of the energy in the 
flow as a function of F at different times. Most of the energy in the flow is within a narrow range in F, of less than a factor 
of 2, around {T}e, mcd- Note that {T)E,inad is also very close to (F)*, which is {r)_B calculated over the region ^ > and is 
close to {F)^; calculated over the whole flow at t < t^ — to^o- It can be seen that at the time when {r)E,med ~ F(5«) (after 
which we can no longer calculate (F)_B.mcd semi-analytically) , we have {F)_B.mcd = F(^, ) ~ ao, and the Lorentz factor of the 
plasma in the region ^ > i^* is between ~ ctq and ~ 2cro. Since at that stage most of the magnetic energy is already converted 
into kinetic energy, as (t(^») ~ 0.1 (see Fig. [3}, then this should be close to the asymptotic value of {F)£;.med at late times. 
Therefore, we can see that all along, from early to late times, {r)_B,incd is very close to {F)_b. This supports the choice of {F)_b 
as being representative of the typical Lorentz factor of the material carrying most of the energy in the lab frame. 

The right panel of Fig. IG2I shows a calculation in the CM frame of the cumulative energy (E') in the flow as a function 
of its velocity (/3') or Lorentz factor (F'), and supports the picture described in § [4] below Eq. (|33[l . The evolution of the 
part of the flow ahead of the secondary rarefaction wave (5' > ^*), which is described by an analytic self-similar solution, 
is followed from the time tg when the original rarefaction wave reaches the wall, and the wall is replaced by vacuum. At 
t' = t'o this region covers the whole flow, and it is easy to see that in the CM frame almost half of the total energy (Eq) 
is carried by material with a negative velocity (/?' < 0). This is expected, since by definition the total momentum vanishes 
in the CM frame (P' = 0), and remains so at later times as well, when a good part of the total energy is at ^' < ^l. In 
the CM frame we have {o-'{t' = tg)) ~ ctq''^, since in the bulk of the fiow t[) corresponds to t ~ toao in the lab frame and 
{ait = toao)) ~ (Jo/^ This means that the flow is still highly magnetized at t'o, and subsequently its front part accelerates 
while its back side decelerates, as the magnetic energy is transformed into kinetic energy. At t'o < t' < t'^ ^ t'oao we have 
{r'){a') ~ E'o/Mo — Eo/TcmMo ~ a^^'^ so that the magnetization drops as {a') ~ a^^^ {t' /t'o)~^^'^ while the typical Lorentz 
factor increases as (F') ~ a^^^ {t' /to)^^^ . At t' = t'^ the combined acceleration at the front and deceleration at the back 
saturates as the magnetic and kinetic energies become comparable, {a') ~ 1, and the typical Lorentz factor approaches its 
asymptotic value, (F') ~ . Similar to the lab frame, at t' > t'^ there is also a coasting phase in the CM frame, where 
(F') ~ af/^ while the magnetization continues to drop as {a') ~ t'o/t' . However, in the CM frame there is comparable mass 
and energy in material with F' ~ moving in the positive and negative s'-directions, so that the total momentum adds 
up to zero. 
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